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PREFACE 


This report covers the work completed on the research project "Analysis 
and Computation of Internal Flow Field in a Scramjet Engine," for the period 
ending November 30, 1985. The work was supported by the NASA Langley 
Research Center (Computational Methods Branch of the High-Speed Aerodynamics 
Division) through research grant NAG-1-423. The grant was monitored by Dr. 
A. Kumar and Mr. J. R. Drummond of the High-Speed Aerodynamics Division. 
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ABSTRACT 


Analyses and numerical procedures are presented to investigate the 
radiative interactions in transient energy transfer processes in gaseous 
systems. The nongray radiative formulations are based on the wide-band 
model correlations for molecular absorption. Various relations for the 
radiative flux are developed; these are useful for different flow conditions 
and physical problems. Specific plans for obtaining extensive results for 
different cases are presented. The methods presented in this study can be 
extended easily to investigate the radiative interactions in realistic flows 
of hydrogen-air species in the scramjet engine. 
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1. INTRODUCTION 


In the past two decades, a tremendous progress has been made In the 
field of radiative energy transfer in nonhomogeneous nongray gaseous systems. 
As a result, several useful books [1-18] and review articles [19-26] have 
become available for engineering, meteorological, and astrophysical 
applications. In the sixties and early seventies, radiative transfer analyses 
were limited to one-dimensional cases. Multidimensional analyses and 
sophisticated numerical procedures emerged in the mid-to-late seventies. 
Today, the field of radiative energy transfer in gaseous systems is getting an 
ever increasing attention because of its applications in the areas of the 
earth's radiation budget studies and climate modeling, fire and combustion 
research, entry and reentry phenomena, hypersonic propulsion and defense- 
oriented research. 

In most studies involving combined mass, momentum, and energy transfer, 
the radiative transfer formulation has been coupled only with the steady 
processes. The goal of this research is to include the nongray radiative 
formulation in the general unsteady governing equations and provide the step- 
by-step analysis and solution procedure for several realistic problems. The 
specific objective of the present study is to investigate the one-dimensional 
transient radiative transfer in a nongray gaseous system. In the future work, 
the present analysis will be extended (in a systematic manner) to the problems 
of combined transfer processes in chemically reacting flows. 

For the present study, the information on band absorption and correlation 
is summarized in section 2 and fundamental radiative flux equations are 
presented in section 3. The basic formulation for the transient radiation is 
given in section 4, and this is applied to a special case in section 5. The 
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solution procedures are described In section 6, and plans for obtaining 
specific results are presented In section 7. 

2. BAND ABSORPTION AND CORRELATIONS 

The study of radiative transmission In nonhomgeneous gaseous systems 
requires a detailed knowledge of the absorption, emission, and scattering 
characteristics of the specific species under Investigation. In absorbing and 
emitting mediums, an accurate model for the spectral absorption coefficient Is 
of vital Importance In the correct formulation of the radiative flux 

equations. A systematic representation of the absorption by a gas. In the 

Infrared, requires the Indentlflcatlon of the major infrared bands and 
evaluation of the line parameters (Tine Intensity, line half-width, and 
spacing between the lines) of these bands. The line parameters depend upon 
the temperature, pressure and concentration of the absorbing molecules and. In 
general, these quantities vary continuously along a nonhomogeneous path In the 
medium. In recent years, considerable efforts have been expended In obtaining 
the line parameters and absorption coeffi cents of Important atomic and 

molecular species [27-30]. 

For an accurate evaluation of the transmittance (or absorptance) of a 
molecular band, a convenient line model Is used to represent the variation of 
the spectral absorption coefficient. The line models usually employed are 
Lorentz, Doppler, and Voigt line profiles. A complete formulation (and 
comparison) of the transmittance and absorptance by these line profiles Is 
given In [22-26]. In a particular band consisting of many lines, the 

absorption coefficient varies very rapidly with the frequency. Thus, It 
becomes very difficult and time-consuming task to evaluate the total band 
absorptance over the actual band contour by employing an appropriate line 
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profile model. Consequently, several approximate band models (narrow as well 
as wide) have been proposed which represent absorption from an actual band 
with reasonable accuracy [22-26, 31-40]. Several continuous correlations for 
the total band absorption are available in literature [22-26, 36-40], These 
have been employed in many nongray radiative transfer analyses with varying 
degree of success [22-26, 41]. A brief discussion is presented here on the 
total band absorption, band models, and band absorptance correlations. 

The absorption within a narrow spectral interval of a vibration rotation 
band can quite accurately be represented by the so-called "narrow band 
models." For a homogeneous path, the total absorptance of a narrow band is 

given by 

- / [l-exp(k X)] dw (2.1) 

Ao) 

where is the volumetric absorption coefficient, u is the wave number, 

and X = py is the pressure path length. The limits of integration in Eq. 
(2.1) are over the narrow band pass considered. The total band absorptance of 
the so-called "wide band models" is given by 

00 

A = / [l-exp(-k^X)] d(a)-to^) (2.2) 

««oo 

where the limits of integration are over the entire band pass and is the 
wave number at the center of the wide band. In actual radiative transfer 

analyses, the quantity of frequent interest is the derivative of Eqs. (2.1) 
and (2.2). 

Four commonly used narrow band models are Elsasser, Statistical, Random 
Elsasser, and Quasi -Random. The application of a model to a particular case 
depends upon the nature of the absorbing emitting molecule. Complete 
discussion on narrow bands models, and expressions for transmittance and 
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integrated absorptance are available in the literature [22-26, 31-33]. 

Detailed discussions on the wide band models are given in [22-26, 34-40]. The 
relations for total band absorptance of a wide band are obtained from the 
absorptance formulations of narrow band models by employing the relations for 
the variation of line intensity as [22-26, 37-40] 

Syd = (S/AQ)exp{[-bQ|w-a)Qi]/AQ} {2.3) 

where Sj is the intensity of the jth spectral line, d is the line spacing, S 
is the integrated intensity of a wide band, Ag is the band width parameter, 
and bg = 2 for a symmetrical band and bg = 1 for bands with upper and lower 
wave number heads at The total absorptance of an exponential wide band, 

in turn, may be expressed by 

A(u,P) H A(u,P)/A^ = ^ [Ayu,p)]d(a)-a)Q) (2.4) 

° band 

where u * SX/Ag is the nondimensional path length, p = is the line 

structure parameter, Y|_ is the Lorentz line half-width, and Ajj(u,P) repre- 
sents the mean absorptance of a narrow band. 

By employing the Elsasser narrow band absorptance relation and Eq. (2.3) 
the expression for the exponential wide band absorptance is obtained as 
[25,16] 

% 

A(u,p) = Y + (l/ii) / [In 4> + E^ Ml dz (2.5) 

0 

where 4> * u sinh p/(cosh p- cos z), y “ 0.5772156 is the Euler's constant, 
and Ej^{(|>) is the exponential integral of the first order. Analytic solution 

of Eq. (2.5) can be obtained in a series form as [25, 26] 

00 

A(u,p) = S {-{A)"[SUM(mn)]/[n(B+l)"n!(n-l)!]} 
n=l 


( 2 . 6 ) 
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where 

SUM (mn) = S C(n+m-l) ! (2m-l) !C"’]/(2"’(ml )^] 
m=o 

A = - u tanhp, B = l/coshp, 

C = 2/(l+coshp) - 2B/(B+1). 

The series in Eq. (2.6) converges rapidly. When the weak line approximation 

for the Elsasser model is valid (i.e. p is large), then Eq. (2.5) reduces to 

A(u) = Y + in(u) + Ej^(u). (2.7) 

In the linear limit, Eqs. (2.5) and 12.6) reduce to A * u, and in the 

logarithmic limit they reduce to A - y + ln(u). It can be shown that Eq. 
(2.5) reduces to the correct limiting form in the square-root limit. Results 
of Eqs. (2.5) and (2.6) are found to be identical for all pressures and 

pathlengths. For p > 1 atm, results of Eqs. (2.5)-(2.7) are in good agreement 
for all path lengths. 

By employing the uniform statistical, general statistical, and random 

Elsasser narrow band models absorptance relations and Eq, (2.3), three 

additional expessions for the exponential wide band absorptance were obtained 
in [25, 26]. The absorptance results of the four wide band models are 

discussed in detail in [26]. The expression obtained by employing the uniform 
statistical model also reduces to the relation (2.7) for large p. 

Several continuous correlations for the total absorptance of a wide band, 
which are valid over different values of path length and line structure 
parameter, are available in the literature. These are discussed, in detail, 
in [22-26, 37-40] and are presented here in the sequence that they became 

available in the literature. Most of these correlations are developed to 
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satisfy at least some of the limiting conditions (nonoverlapping Tine, llenar, 
weak line, and strong line approximations, and square-root, large pressure, 
and large path length limits) for the total band absorptance [23-26]. Some of 
the correlations even have experimental justifications [22,35]. 

The first correlation for the exponential wide band absorptance (a three 
piece correlation) was proposed by Edwards et al. [34, 35]. The first 
continuous correlation was proposed by Tien and Lowder [22], and this Is of 
the form 


Alu,p) = ln(uf(t){u+2)/[u+2f(t)]}+l) (2.8) 


where 


f(t) = 2.94[l-exp(-2.60t)], t = 6/2. 

This correlation does not reduce to the correct limiting form In the square- 
root limit [23,26], and Its use should be trade for p > 0.1. Another 
continuous correlation was proposed by Goody and Belton [39], and In terms of 
the present nomenclature, this Is given by 

A(u,6) = 2 ln{l+u/[4+(itu/4t)]^^^}, p * 2t. (2.9) 

Use of this corelation Is restricted to relatively small p values [23-26]. 
Tien and Ling [40] have proposed a simple two parameter correlation for 
A(u,p) as 

A(u) * s1nh“^ (u) (2.10) 

which Is valid only for the limit of large p. A relatively simple 

continuous correlation was Introduced by Cess and Tlwarl [23], and this Is of 


the form 
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A(u,p) = 2 lna+u/{2+Cu{l+l/p)]^/^}) (2.11) 

where p = 4t/n = 2p/w. By slightly modifying Eq. (2.11), another form of 
the wide band absorptance is obtained as [25, 26] 

A(u,p) = 2 ln(l+u/{2+[u(c+ir/2p)]^^^}] (2.12) 


where 


S O.l, p < 1 and all u values 

0.1, p > 1 and u < 1 

0.25, p > 1 and u > 1. 


Equations (2.11) and (2.12) reduce to all the limiting forms [23]. Based on 
the formulations of slab band absorptance, Edwards and Balakrishnan [37] have 
proposed the correlation 

A(u) = ln(u) + E^(u) + Y +1 - £3(0) (2.13) 

which is valid for large p. For present application, this correlation 

should be modified by using the technique discussed in [25, 26]. Based upon 
the formulation of the total band absorptance from the general statistical 
model, Felske and Tien [38] have proposed a continuous correlation for 
A(u,p) as 

A(u,p) = 2E^(tPy) + ^ 1 ^ 9 J 2) - E^[(Py/2)(l+2t)] 

+ ln[tpy)^/(l+2t)] + 2 y (2.14) 

where 

Py = {(t/u)[l + (t/u)]}' 
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The absorptance relation given by Eq. (2.7) is another siinpTe correlation 
which is valid for all path lengths and for t * (p/2) > 1. The relation of 
Eq. (2.6) can be treated as another correlation applicable to gases whose 
spectral behavior can be described by the Elsasser model. In [26] Tiwari has 
shown that the Elsasser as well as random band model formulations for the 

total band absorptance reduce to Eq. (2.7) for t > 1. 

Band absorptance results of various correlations are compared and 
discussed in some detail in [25, 26* 41]. It was found that results of these 

correlations could be in error by as much as 40% when compared with the exact 

solutions based on different band models. Felske and Tien's correlation was 
found to give the least error when compared with the exact solution based on 
the general statistical mode! while Tien and Lowder's correlation gave the 

least error when compared with the exact solution based on the Elsasser model. 
The results of Cess and Tiwari Vs correlations followed the trend of general 
statistical model. Tiwari and Batki's correlation [Eq. 2.6 or 2.7] was found 
to provide a uniformly better approximation for the total band absorptance at 
relatively high pressures. The sole motivation in presenting the various 
correlations here is to see if their use in actual radiative processes made 
any significant difference in the final results. 

In reference 41, use of several continuous correlations for total band 
absorptance was made to two problems to investigate their influence on the 
final results of actual radiative processes. For the case of radiative 
transfer in a gas with internal heat source, it was found that actual center- 
line temperature results obtained by using the different correlations follow 
the same general trend as the results of total band absorptance by these 
correlations. From these results, it may be concluded that -use of the Tien 
and Lowder's correlation should be avoided at lower pressures, but its use is 
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justified (at moderate and high pressures) to gases whose spectral behavior 
can be described by the regular Elasasser band model. For all pressures and 
path length conditions, use of the Cess and Tiwari's correlations could be 
made to gases with bands of highly overlapping lines. In a more realistic 
problem Involving flow of an absorbing emitting gas, results of different 
correlations (except the Tien and Lowder's correlation) differ from each other 
by less than 6% for all pressures and path lengths. Use of Tien and Lowder's 
correlations Is justified for gases like CO at moderate and high pressures. 
For gases like C02» use of any other correlation Is recommended. While Felske 
and Tien's correlation Is useful for all pressures and path lengths to gases 
having random band structure. Tlwarl and Batkl's simple correlation could be 
employed to gases with regular or random band structure but for P >1.0 atm. 

3. RADIATIVE FLUX EQUATIONS 

For many engineering and astrophysical applications, the radiative 
transfer equations are formulated for one-d1mens10nal planar systems (Fig. 
3.1). For diffuse boundaries and 1n the abscence of scattering, expressions 
for the radiative flux and Its derivative are given as [8] 


+ 2[J ^ e|^j^(t) ^ ■ ■i^^)dt] 


•br 


(3.1) 


and 


^ ^IX ^ ®2X ^2^'^oX ~ '^X^ 


+ 2[/ ■ tildt 


(3.2) 



9a 



Figure 3.1 Plane radiating layer between parallel boundaries. 
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where 

y L 

\ ■ J ‘'y • \\ ■ I 
0 0 

E^{t) = dll 


(3.3a) 

(3.3b) 


In the preceding equations, exponential integral functions, and 

and represent the optical coordinate and optical path, respectively. 
The quantitites and represent the spectral surface radios! ties and 
for nonreflecting surfaces, B^j^ - e^^ = etc. Thus, for non- 

reflecting boundaries, Eqs. (3.1) and (3.2) are expressed in terms of the wave 
number as (see Appendix A) 

‘ ®lu ‘ ®2o) 


T T 

+ 2 C/ “ Pyt) ^2^^^ h (3.4) 

0 w 

and 

T T# 

+ 2 [/ “ F, (t) E, (T^-t)dt + F2^(t) E.(t-ydt3 (3.5) 

0 0) 

where 




Equations (3.4) and (3.5) are the general equations for one-dimensional 
absorbing-emitting medium with diffuse non-reflecting boundaries. For nongray 
analyses, it is often convenient to replace the exponential integrals by 
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appropriate exponential functions [6, 8]. Upon employing the exponential 
kernal approximation [8] 




Eqs. (3.4) and (3.5) are expressed in physical coordinatates as 
%0) (y) = ®lco - ®2a, 

I ^0 '^0) I '^0) 

- I I ^20) *^0) I ^3.6) 

0 

■" \ ''20)^^^ I '^0) (Z“y)3dz (3.7) 

where z is a dummy variable for y. However, by differentiating Eq. (3.6) 
directly, there is obtained 

''la "a I 

0 

I ^y '"2(0 ''w ®''P'^' I >«u,{2-y)]dz 


(3.8) 
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The slight difference in Eqs. (3.7) and (3.8) should be noted. This is a 
consequence of using the exponential kernal approximation. If one has to make 
a decision as which equation to use* it is recommended to use Eq. (3.8). 

The total band absorptance, as given by Eq. (2.2), can be expressed in a 
slightly different form as 

” -1 

A(y) = / [1 - exp(- k^ y)] du ~ cm" (3.9a) 

0 

where both k^ and u have units of cm“^. Differentiation of Eq. (3.9a) gives 

“ 

A' (y) * J k^ exp(- y) du ~ cm (3.9b) 

0 

and 

CO 

A"(y) = / - k^ expl- kj^ y) do) ~ cm”^ (3.9c) 

0 

Equations (3.9) are employed to express Eqs. (3.6) and (3.8) in terms of the 
band absorptance. 

The total radiative flux is given by 

00 

QR{y) = / (3.10) 


such that 


dq,^(y) 

dy 


r '“'ru 


du) 




(3.11) 
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Upon substituting Eq. (3.6) into Eq. (3.10) and Eq. (3.8) into Eq. (3.11) 
there is obtained for a multiband gaseous system 

qR(y) = ej - ej 


1=1 Awi 0 

L . 

* ^ '^ui 7 *^u-i^2-y)]dz}da). (3.12) 


dqf^(y) 

3y 


1=1 A(i)l 


I . exp[-|k„,(y-z)]dz 

1 = 1 Ao)l 0 


+ /y '^2o)i^^^ I ka^i(z-y)3dz}dw. 


(3.13) 


It should be pointed out that the following relations have been used in 
obtaining Eqs. (3.12) and (3.13) 


I e^^ do) = ep / eg^ da> = eg 
0 0 


/ {/ expC- I kJy-z)]dz}do) 

0 0 


" t . ''.oi 7 '‘^>'^“1 

1 = 1 Awi 0 


where n represents the number of bands in a multiband system. 
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By utlizing the deffnttfons of the band absorptance and its derivatives 
as given by Eqs. (3.9) and evaluating the value of the Planck function at the 
center of each band, Eqs. (3.12) and (3.13) are expressed as 

qR (y) » ei - eg 


n ,y 


i*l 0 ^"oi ’ 2 


■ ^ a;. Ci (z-y)]dz} 

y -oi ^ ^ 


dqR(y) , n 

1=1 01 01 AcJl 


(3.14) 




^ ^lo) A'i 4 (y-z)idz 
^ i=l o %i ^ ^ 


+ / (z) A*;, cl (z-y)]dz} 

y ‘^“oi ^ ^ 


(3.15) 


where represents the center of the ith band. 

Equations (3.14) and (3.15) are in proper form for obtaining the nongray 
solutions of molecular species. However, in order to be able to use the band 
model correlations, these equations must be transformed in terms of the 
correlation quantities defined in' Eq. (2.4). The following quantities, 
therefore, are needed for the transformation 

u = (S/Aq) py; Uq = (S/A^) PL; PS = / do) (3.16) 

Ao) 


Now, by using the definition A = A/A^, Eq. (3.9b) is written as 
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A' (y) = ^ [A(y)] = ^ = A' (y)/A^ ~ cm'^ 

Thus, 


A'(y)=A„i^.A„C^^] = PS(T) A'(u) 


(3.17a) 


Similarly 

A“(y) = [P S(T)]2 (1/Ajj) A"(u) (3.17b) 

The dimensions of both sides in Eqs. (3.17a) and (3.17b) agree with the 
dimensions given in Eqs. (3.9b) and (3.9c). By employing the definitions of 
Eqs. (3.16) and (3.17), Eqs. (3.14) and (3.15) are expressed as 

q„(u) • ej - ej 


dqj^(u) 

du 


u. 


^ 7 4 ^oi '^10). ^ ^"1 


0 1 


- 1 f2». <“i> *i '“i - '‘"P 


>\ »• r3 /..I 


Ui - 1 


I *01 


u- 




0 1 


u . 
01 


+ / 1^20).^*^'^ ^'i 


u. 1 
1 


(3.18) 


(3.19) 


where u' is the dummy variable for u and A'(u) = dA/du. It should be noted 

that and in Eqs. (3.18) and (3.19) represent the values of 

i i 
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and Fg^ at the center of the ith band, and dqj^/dy = (dqj^/du){du/dy) = 
[P S(T)/A^] (dq,^/du). 

By defining the new independent variables as 

I - u/Uq = y/L; V = u7Uq = z/L (3.20) 

Eqs. (3.18) and (3.19) can be expressed as 


“ ®1 ” ®2 


*1.", *01 %i 

1=1 0 1 




j " ^<5)5 <*oi %1» 


n ^ 

+ 1 i A„, Um </ F„„ (5') A'J [|o„,(5-r)]d5' 


T "01, “oi 


*i ^20.1 («■) S'S C|o„,(5’-5)]d5’) 


(3.21) 


(3.22) 


where again A' (u) dentoes the derivative of A(u) with respect to u, and 
dq,^/du = (dqj^/d5) (d?/du) = (1/u^) (dq,^/dC). 

Equations (3.18) through (3.22) allow us to make use of the band model 
correlations for the wide-band absorptance because these correlations are 
expressed in terms of u and p. However, it is often desirable and 
convenient to express the relations for qj^ and div q[^ which only involve 
A(u) and A'(u) but not A''(u). This is accomplished by integrating the 
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integrals in the expressions for qj^ and div by parts. This results in 
simpler integrals. Upon performing the integration by parts on the integrals 
in Eqs. (3.21) and (3.22), the equations can be expressed in alternate forms 
as (see Appendix B) 

q„(5) - ej - 


n 

+ 2 
i=l 


dq,^ (C) _ 3 

“ 1 




* / 


n 

t A . 


1»1 


01 


Cde,, (O/d?'] A. C| u^.(5-C') dV 

0)^ 1 4 01 

[de^ (C')/dC3 A. C|u^. iVOldV} 
u . {/ [de (C' )/dC'] A', \x Al-V )1dl' 

U1 ^ u)^ 1 4 01 

- / Cde^ (C')/dC'] a;. [|uQ.(C'-^)]dC'} 
C i 


(3.23) 


(3.24) 


It should be noted that Eq. (3.24) can be obtained directly by differentiating 
Eq. (3.23) with respect to C using the Leibnitz formula. This is shown in 
Appendix B. 

Equations(3.21) , (3.23), and (3.24) are the most convenient equations to 
use when employing the band-model correlations in radiative transfer analyses. 
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4. BASIC FORMULATION FOR TRANSIENT PROCESSES 

The interaction of radiation in transient transfer processes has received 
very little attention in the literature. Yet, the transient approach appears 
to be the logical way of formulating a problem in general sense for- elegant 
numerical and computational solutions. The steady-state solutions can be 
obtained as limiting solutions for Targe times, 

A few studies available on radiative ineractions reveal that the 
transient behavior of a physical system can be influenced significantly in the 
presence of radiation [42-45], Lick investigated the transient energy 
transfer by radiation and conduction through a semi-finite medium [42], A 
kernal substitution technique was used to obtain analytic solutions and 
display the main features and parameters of the poblem, Doornink and Hering 
studied the transient radiative transfer in a stationary plane layer of a 
nonconducting medium bounded by black walls [43]. A rectangular Milne- 
Eddington type relation was used to describe the frequency dependence of the 
absorption coefficient. It was found that the cooling of the layer initially 
at a uniform temperature is strongly dependent on the absorption coefficient 
model employed, Larson and Viskanta investigated the problem of transient 

combined laminar free convection and radiation in a rectangular enclosure 
[44]. It was demonstrated that the radiation dominates the heat transfer in 
the enclosure and alters the convective flow patterns significantly. The 
transient heat exchange between a radiation plate and a high-temperature gas 
flow was investigated by Melnikov and Sukhovich [45], Only the radiative 
interaction from the plate was considered; the gas was treated as a non- 
participating medium. It was proved that the surface temperature is a 
function of time and of longitudinal coordinate. 
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The objective of this study Is to Investigate the Interaction of nongray 
radiation In transient transfer processes In a general sense. Attention, 
however, will be directed first to a simple problem of the transient radiative 
exchange between two parallel plates. In subsequent studies, the present 
analysis and numerical techniques will be extended to Include the flow of 
homogeneous, nonhomogeneous, and chemically reacting species In one- and 
mul t1 -dimen si onal sys terns . 

The physical model considered for the present study Is the transient 
energy transfer by radiation in absorbing-emitting gases bounded by two 
parallel gray plates (Fig. 4.1). In general, Tj^ and T 2 can be a function of 
time and position and there may exist an Initial temperature distribution in 
the gas. It is assumed that the radiative energy transfer In the axial 
direction is negligible In comparison to that In the normal direction. 

For radiation participating medium, the equations expressing conservation 
of mass and momentum remain unaltered, while the conservation of energy. In 
general. Is expressed as [8] 

P Cp ^ div (k grad T) + PT -^ + U (t> - div q^^ (4.1) 

where p is the coefficient of thermal expansion of the fluid and (j) Is the 
Rayleigh dissipation function. For a semi-infinite medium capable of 
transferring energy only by radiation and conduction, Eq. (4.1) reduces to 


p c 


p at " ay 


(4.2) 


where q Is the sum of the conductive heat flux q = - k (aT/ay) and the 
radiative flux q|^. For the physical model where radiation Is the only mode of 
energy transfer, the energy equation can be written as 
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5T _ 

^ at ~ dy 


(4.3) 


Use of this simplified equation Is made to Investigate the transient behavior 
of a radiation participating medium. 


As pointed out In the previous section (Sec. 3), Eqs. (3.21) and (3.23) 
are convenient equations for the radiative flux. Equations (3.22) and (3.24) 
are two expressslons for the div qR(y), but Eq. (3.24) Is preferred because It 
only Involves the first derivative of A and avoids singularities in the 
large path length limit. 


Upon defining nondlmenslonal radiative heat flux by 


Q (4,t) * q(^(?,t)/[e^(t) - e2 (t)] 


(4.4) 


Eq. (3.21) can be written as 

3 " ^ - 3 

Q il.t) = 1 + - E A . u . {/ c,. (c'.t) a; [i u . (C-C')]dC' 
2 , , 01 01 li 1 2 01 

1=1 0 


- / C-. (C'.t) A! C| u . (5'-C)]d5'} (4.5) 

^ 1 2 01 


where 


(?.t) = F^. (C.t)/[e^(t) - e 2 (t)] 

Equation (4.5) provides the general expression for the radiative flux In the 
nondlmenslonal form. A similar nondlmenslonal form can be obtained also from 
Eq. (3.23). 

By defining ^ (C.t) = T(C,t)/T^ with representing some constant 
reference temperature, Eqs. (4.3) and (3.24) can be combined to yield the 
energy equation In nondlmenslonal form as 
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- a<(> (C,t)/dt =1 z: {/ <1;^^. (5*.t) a; [| {c-5')d5' 

r=l 0 

1 3 

- / 4, {^',t) a;. Ku . (C'-5)]dC*} (4.6) 

- M 1 c 01 

c 

where 

= {PS.(T)[d e^^.{^,t)/55]/(p Cp yt^)} 

The time t tn Eq. (4,6) is defined as t* = t/t^ with representing some 
characteristic time scale of the physical problem; however, for the sake of 
convenience, the asterisk is left out here as well as in further developments. 
From the definitions of 4>(^.t) and (C,t), it should be noted that Eq. 
(4.6) is a nonlinear equation in T(C»t). Equation (4.6), therefore, 
represents a general case of the transient energy by radiation between two 
semi-infinite parallel plates. A similar expression can be obtained also by 
combining Eqs. (4.3) and (3.22). 

5. A SPECIAL CASE OF TRANSIENT INTERACTION 

As a special case, it is assumed that the entire system initially is at 
the fixed (reference) temperature T^. For all time, the temperature of the 
upper plate is maintained at the constant temperature equal to the reference 
temperature, i.e., T 2 = Tq. The temperature of the lower plate is suddenly 
decreased to a lower but constant temperture, i.e., T^ < T 2 . The problem, 
therefore, is to investigate the transient cooling rate of the gas for a step 
change in temperature of the lower plate. 

Since small temperature differences have been assumed and the absorption 
coefficient has been taken as independent of temperature one may comploy 
additionally the linearization. 
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(T-T„) (5.X) 

W 

where again the subscript i refers to the ith band such that is the wave 
number location of the band and represents the temperature of the reference 
wall which could be either Tj^ or T 2 . For the special case considered, since 
we are interested in investigating the transient behavior of the gas because 
of a step change in temperature of the lower plate, is taken to be equal to 
Tj^. Thus, 

e^^i(t,t) - e^^^.{0,t) = (d e^./dT)^^ (T-T^ (5.2a) 


•wi 


(l,t) -e^.(0,t) * (d e^./dT)^^ (Tg-T^) 


(5.2b) 


e^^i(C,t) - e^.(l,t) = (d e^./dT)^ (T-Tg) (5.2c) 

Note that Eq. (5.2c) is obtained by subtracting Eq, (5.2b) from Eq. (5.2a). 
Also, for linearized radiation, 

* 4 T - 3 ij (5.3) 

Thus, Oj = * ®2 * T 2 = ^2 ' ^ 

e^ - eg = 4 a T^ (T^- Tg) . 

It should be pointed out that for a single-band gas, the linearization is 
not required because the temperature distribution can be obtained either by 
combining Eqs. (3.22) and (4.3) or from Eq. (4.6) and the radiative heat flux 
can be calculated from Eqs. (3.21), (3.23), or (4.5). However, for the case 
of multiband gases and for systems involving mixtures of gases, it is 
convenient to employ the linearization procedure in order to use the 
information on band model correlations. The following definition are useful 
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in expressing the governing equations in linearized forms: 

0 = (T-Tji/CTg-Tj.) 


(5.4a) 


•^li ' '^li * '^li = 

“i * V" '=p> '<1 ■ ' jj "ii <=•''=» 

”l1 = ‘V- P 'p> "ll • «11 ' *01 <" 

n 

”l “ ^ V*- P Pp^ ”l » ”l “ ”li 


«ii = “li 


01 


«H = ^11 


(5.4f) 


where H|, Kj^, Nj and represent the values of H, K, N and M evaluated at the 
temperature Ti. As explained in Refs. 8 and 23, these quantities represent 
the properties of the gaseous medium. 

By employing the definitions of Eqs. (5.2) - (5.4), relations for the 
radiative flux, as given by Eqs. (3.21) and (3.23), are expressed as 

Q(C.t) = 1- (3/8 a tJ) Z Uq^. 0(^',t) Al[| u^^. (C-5' ) ]d?‘ 

i *1 0 


1 




5 


Cl - e(r .t)] AJ[| UQ. (5'-C)]d?'} 


(5.5a) 


Q(C.t) = 1 - (1/4 0 tJ) Z A^.[| u^j.{C-5')]dC' 

i=l 0 


and 
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ae(c',t) 

\ ' W 


A.[| u^.(r-5)]dC'} 


(5.5b) 


Thus, the expressions for the heat flux at the lower wall 


Q(O.t) = 1 - {3/8 aT|) Z / Cl-e(?' ,t)]A.(| u^. 

1=1 0 


and 


Q(0 


,t) . 1 . (1/4 c if) Hj, / A,(| 


are given by 

V)dV 

V)dV 


(5.6a) 


(5.6b) 


It should be pointed out that Eqs. (5.5a) and (5.6a) are convenient forms for 
the optically thin and general solutions while Eqs. (5.5b) and (5.6b) are 
useful for solutions In the large path length limit. Once the solutions for 
9(C,t) are known from the energy equation, the appropriate relations for the 
heat flux can be obtained from Eqs. (5.5) and {5.6). 

By employing the definitions of Eqs. (5.2) - (5.4), a combination of 
Eqs. (3.22) and (4.3) provides one form of the energy equation and Eq. (4.6) 
is transformed to obtain another form; these are expressed as 

+ 3 e(e,t) - I Nj = 


? “^ll %i A|'C| UQ.(C-C')]dC’ 

1=1 0 


and 


+ Ce(5',t) - 1 ] A'[[| u^.(C'-C)M'} 


ae(c,t) _ 3 ; ,, ae(?',t) 7.^3 

— aT — 9 ^ — A [-u (c-C')M 

ot 2 . , li oC 1 2 01 

1=1 o 


{5.7a) 


o 



(5,7b) 


- / 




The Initial and boundary conditions for Eq. (5.7) are specified as 

e(C,0) = 1 ; 0(0, t) = 0 ; 0(1, t) = 1 (5.8) 

The parameters in Eq. (5.7) are Nj^ and Uq. For a given gas, the parameters 
are the gas pressure and the temperature of the lower wall. Equation (5.7b) 
is the convenient form for solutions in the Targe path length limit. 

6. METHOD OF SOLUTIONS 

For the general case, the temperature distribution is obtained from the 
solution of the energy equation, Eqs. (5.7). Once 9(5, t) is known, the 
radiative heat flux is calculated by using the appropriate form of Eq. (5.6). 
Before discussing the solution procedure for the general case, however, it is 
desirable to obtain the limiting forms of Eqs. (5.5) and (5.7) in the 
optically thin and large path length limits and investigate the solutions of 
resulting equations. 


6.1 Optically Thin Limit 

In the optically thin limit [8, 23], A(u), = u, AV(u) * 1, and 

A"(u) = 0. In this limit, therefore, Eq. (5.7a) reduces to 

+ 3 Nj 0(5,t) - I Nj = 0 (6.1a) 

From an examination of Eq. (6.1a) along with the definitions given in Eq. 
(5.4), it is evidient that in the optically thin limit the temperature distri- 
bution in the medium is independent of the 5- coordinate. This is a 
charateristic of the optically thin radiation [8, 23]. Thus, Eq. (6.1a) can 
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be written as 

3 N,(t) e(t) - |n ft) = 0; e(C.O) = 1 (6.1b) 

dt 1 2 1 

Since gas properties are evaluated at known reference conditions, is 
essentially constant, and solution of Eq. (6.1b) is found to be 

9(t) = |[1 + exp(- 3 Njt)] (6.2) 


In the optically thin limit, Eq. (5.7b) reduces to 


S9(?,t) . 3 , " „ ^ mV.t) 


5t 


(6.3a) 


0 ' I 

A differentiation of Eq. (6.3a) with respect to C (by using the Leibnitz's 
rule) results in 


d 


+ 3 0(c,t)] = 0 


(6.3b) 


or 

+ 3 Nj e(5,t) = C(t) (6.3c) 

The constant of integration C(t) is evaluated through the combination of Eqs. 
(5.8) and (6.3a) and is found to be C(t) “ A substitution of this in 
Eq. (6.3c) gives Eq. (6.1a) for which the solution is given by Eq. (6.2). 
Thus, as would be expected, both general forms of the energy equation reduce 
to the same equation in the optically thin limit. 

In the optically thin limit Eqs. (5.5a) and (5.5b) respectively reduce to 
, 5 1 

Q(?,t) = 1 - C3/(8aTf)] (PLK.) {/ 9(?',t)dC' + J Cl-9(C',t)]dC'} (6.4a) 

^0 I 


and 
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QU.t) = 1 - [3/8orTp3 (PLK^) {/ 


d9 

?>V 


(^-^')dC’ + / 


99 

W 


(C'-e)d^'} (6.4b) 


Through integration by parts, it can be shown that Eq. (6.4b) reproduces Eq. 
(6.4a). By noting that, in the optically thin limit, 9U,t) = 9(t), Eqs. 

(6.4) can be expressed as 

Q(C,t) * 1 - C3/(8cjl|)] (PLKj) C(l-C) + (2C-1) 9(t)3 (6.5) 

It should be pointed out that Eq. (6,5) can be obtained directly from Eq. 
(6.4b) without performing the inegration by parts. The heat transfer from the 
lower surface in the optically thin limit, therefore, is given by 

Q(0,t) = 1 - C3/(8 cjT^)] (PLKj)[l - 9(t)] (6.6) 

The result of Eq. (6.6) can be obtained directly by letting C “ 0 in either 
of Eqs. (6.4). The relation for 9(t) in Eq. (6.6) is obtained from Eq. 
(6.2). Thus, evaluation of the temperature distribution and radiative heat 

flux in the optically thin limit does not require numerical solutions. 


6.2 Large Path Length Limit 


In the large path length limit (i.e., for Uq^ >> 1 for each band), one 
has A(u) = An(u), A'(u) = 1/u, and A''(u) * - 1/u^ [8, 23]. Thus, in the 

large path length limit, Eq. (5,7a) reduces to 


59(C,t) 
“Tt 


+ 3 9U,t) - I Nj 


I 

= {/ 9(5', t) 


dg' 

(g-g')^ 



Ce(5',t) - 1 ] 


dg’ ^ 
— — 

(g'-g)^ 


(6.7a) 


It should be noted that for any fixed value Tj^ and a given gas, Nj^ and are 
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constants; but, 9(C»t) does depend on 1 . For a given gas and with known 
values of T 2 and tj„, the solution of Eq, (6.7a) can be obtained by specifying 
Tj^. Equation (6.7a) involves singular integrals with Cauchy type kernals and, 
therefore, a closed form solution does not appear to be possible; numerical 
solutions, however, can be obtained by the variation of parameter technique. 
Because of the singular nature of integrals, Eq. (6,7a) is not a convenient 
equation for the large path length limit solutions. 


99(C,t) _ „ 56(5'. t) dV 

— — - “^1 J ■ — ^ — TFTT 


(6.7b) 


In the large path length limit, Eq. (5.7b) reduces to 

1 

f 

0 

Equation (6.7b) is a convenient form for solution in the large path length 
limit. An analytical solution of Eq, (6.7b) may be possible, but numerical 
solution can be obtained quite easily. 


In the large path length limit, Eqs, (5,5a) and (5.5b) reduce 
respectively to 


Q(5.t) = 1 - (1/4T^) H^[/ 9(5', t) 

0 


dV / d?' ^ 

IC-C) ■ I?-?')-' 


(6.8a) 


and 


Q(C,t) = 1 - (l/ 4 alf) S H {/ ^ lyH An d u (C-^')]d^' 

i , , n c 01 

r=l 0 


^ S9(5',t) 

+ 5 ^ 


An[| Ug^ CC'-OldC'} 


(6.8b) 


The expressions for dimensionless radiative heat flux from or to the wall 
are obtained by setting C = 0 in Eqs. (6.8) as 
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Q{0,t) = 1 - (Hj/40t[) 


/ Cl - e(c',t)] #' 
0 ^ 


(6.9a) 


and 

n 1 

Q(0,t) = 1 (l/4aTj) / M|^ An(| u^. V )dV (6.9b) 

Thus, once the temperature distribution is known from solutuions of Eq. (6.7), 
the wall heat flux can be calculated by using the corresponding form of Eqs. 
(6.9). 


6.3 Numerical Solutions of Governing Equations 

General solutions of Eqs. (5.7a) and (5.7b) are obtained numerically by 
employing the n«thod of variation of parameters. For this, a polynomial form 
for 9(C,t) is assumed in powers of t with time dependent coefficients as 


e(C.t) 


n 

S 

m=o 


c (t) 
m 


( 6 . 10 ) 


By considering only the quadratic solution in C, and satisfying the boundary 
conditions of Eq. (5.8), one finds 

e(C,t) = 5^ + g(t) (6.11) 


where g(t) represents the time dependent coefficient. At t = 0, a combination 
of Eqs. (5.8) and (6.11) yields the result 

g(0) = (l-C^)/(?-C^) (6.12) 


Also, from Eq. (6.11) there is obtained 


dt 


dt 


(6.13) 
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and 

= 2 C + g(t) (1-2?) (6.14) 

0 ? 

Equations (6.11) - (6.14) are employed to obtain specific solutions of Eqs. 
(5.5) and (5.7). 

By substituting Eqs. (6.11) and (6.13) in Eq. (5.7a), there is obtained 

g'(t) + G^(?) g(t) = G2(?) (6.15) 

where the integral functions G|(?) and G2(?) are defined in Appendix C. The 
solution of Eq. (6.15) is given by 

g(t) = c expC- Gj(?)ti + G2(?)/Gj(?) (6.16a) 

Since at t = 0, g(t) - g(0), then c = g(0) - G2(?)/Gj^(?). Thus, Eq. (6.16a) 
becomes 

g(t) = [g(0) - G2(?)/G^(?)i exp[- Gi(?)t] + G2(?)/G^(?) (6.16b) 

where g(0) is given by Eq, (6.12). The integrals in functions G^(?) and 
G2(?) can be evaluated easily by numerical means, after substituting the 
relation for A' ' (u) . 

A substitution of Eqs. (6,11), (6.13) and (6.14) into Eq. (5.7b) results 
in 

g'(t) + Ggi?) g(t) = G4(?) (6.17) 

where the integral functions S3(?) and G^(?) are defined in Appendix C. The 
solution of Eq. (6.17) is found to be 

g(t) = Cg(0) - G^ (?)/G3(?)] expC- G3(?)t] + G4(?)/G3(?) 


(6.18) 
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where again g(0) is given by Eq. (6.12). 


The solutions of Eqs. (5.7a) and (5.7b) can be expressed in a convenient 
form as 


e(c,t) 




(6.19) 


In Eq. (6.19), gi(t) is given by Eq. (6,16b) and is used for the solution of 
Eq. (5.7a) and g2(t) is given by Eq. (6.18) and is used in obtaining the 
solution of Eq. (5.7b). The both approach should result in the same final 
solution. 

For the steady state case, the solution again is given by Eq. (6.19), but 
functions gi(t) and g2(t) are no longer a function of time and are given by 

g^ = 62(0/61(5); 92 - 24(0/63(5) (6.20) 

The solutions for the steady case are available in the literature and are 
useful in comparing the results of this study in the limit of t ■* % 

The expressions for the nondimensional radiative flux are obtained from a 
combination of Eqs. (5.5), (6.11) and (6.14) such that 

Q(C.t) = 1 - 65(5) gi{t) - 6g(5) (6.21a) 

and 

Q(C,t) = 1 - 67(5) g2(t) - 6g(5) (6.21b) 

where 6g(C) through 6g(0 are defined in Appendix C, and qi(t) and q2(t) 
are given respectively by Eqs. (6.16b) and 6.18). Consequently, the 
expressions for the radiative heat flux at the lower wall are obtained as 
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Q(0,t) = 1 + Gg q^Ct) - Gjq (6.22a) 

and 

Q(0,t) = 1 + ggCt) + G^2 (6.22b) 

where Gg through Gj^g defined in Appendix C and are not function of C. 
It should be noted that the solutions presented in Eqs, (6,21) and (6.22) 
require the solution of the energy equation as given by Eqs. (6.19) 

6.4 Numerical Solutions of Large Path Length Equations 

As mentioned earlier, Eqs. (6.7b), (6.8b) and (6.9b) are the most appro- 
priate equations to use in the large path length limit. However, numerical 

procedure is presented for both forms of the energy and radiative flux 

equations. Once again Eqs. (6.11) through (6.14) provide the basis for 

numerical solutions also in the large path length limit. For this limit, the 

solution given by Eq. (6.19) is expressed as 


e(C.t) = ^ + [gjjt)’] 


(6.23) 


where g 3 (t) is used for the solution of Eq. (6.7a) and g 4 (t) for Eq. (6.7b). 

A substitution of Eq. (6.23) into Eq, (6,7a) results in 

gj(t) + Gj 3 (C) gjit) = G^4(5) (6.24) 

where integral functions Gj^ 3 ( 5 ) and Gj^^(C) are defined in Appendix C. The 
solution of Eq. (6.24) is found to be 

g 3 (t) = Cg(0) - G^4(C)/G^3(5)] expC- Gj 3 (C)t] + G^4(C)/Gj3(5) 


(6.25) 
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where g(0) is defined again by Eq. (6.12). Equation (6.23) along with Eq. 
(6.25) provides the solution of the energy equation, Eq. (6.7a). 

A combination of Eqs. (6.23) and (6.7b) results in 

9^ ( 94^^) * ^16^^^ (6.26) 

The integral functins 6^5(0 and Gj^gl?) appearing in Eq. (6.26) are defined 
in Appendix C. These, however, can be evaluated easily with the results 

615(5) = CMj/(5-5^)i{2 + (25-1) AnC(5-l)/5)} (6.“27a) 

Gjg(5) * 2 CMj/(5-5^)]{l + 5 An[(5-l)/5]} (6.27b) 

The solution of Eq. (6,26) is found to be 

g4(t) = Cg(0) - Gjg(5)/Gj5(5)] expC- Gj5(5)ti + G^g(l)/G^^(a (6.28) 

where again g(0) is defined by Eq. (6.12). A combination of Eqs. (6.23), 
(6.27) and (6.28) provides the solution of the energy equation, Eq. (6.7b). 
The only parameter appearing in the solution of Eq. (6.7b) is M^. 

The expessions for the nondimensional heat flux in this case is obtained 
from a combination of Eqs. (6,8) and (6.23) as 

Q(5,t) = 1 - Gjy(5) g3(t) - G^g(5) (6.29a) 

and Q(5,t) = 1 - G^g(5) g4(t) - G2 q( 5) (6.29b) 

where Gj^^(5) through GgQ(5) are defined in Appendix C and can be evaluated 
in closed forms. The corresponding expressions for the radiative heat flux at 
the lower wall are found from Eqs, (6.9) and (6.23) as 


Q(0,t) = 1 - G21 g3(t) - G 


(6.30a) 
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and Q( 0 ,t) = 1 - G23 g^Ct) - G24 { 6 . 30 b) 

where again G21 through G24 are defined in Appendix C and are not functions of 
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7. RADIATIVE INTERACTION IN LAMINAR FLOWS 
The physical systen considered is the energy transfer in laminar, in- 
compressible, constant properties, fully-developed flow of absorbing-emit- 
ting gases between parallel plates (Fig, 7,1). The condition of uniform 
surface heat flux for each plate is assumed such that the temperature of the 
plates varies in the axial direction. Fully developed heat transfer is 
considered, and axial conduction and radiation is assumed to be negligible 
as compared with the normal components. Consistent with the constant prop- 
erties flow, the absorption coefficient is taken to be independent of temp- 
erature and radiation can be linearized. Extensive treatment of this prob- 
lem is available in the literature [23, 41]. The primary motivation of 
studying the problem here is to investigate the extent of radiative inter- 
action for high temperature flow conditions. 

7, 1 Basi c Formu 1 at i on 

For the physical problem considered, the energy equation, Eq. (4.1), 
can be expressed as [8J 

pC (II + u H + V II) = k l!l + gTu ^ + y{ii3^-divq„ (7.1) 

^ 9t ax ay ay2 dx ay 

where u and v denote x and y components of velocity, respectively. In de- 
riving Eq. (7.1) it has been assumed that the net conduction heat transfer 
in the x direction is negligible compared with the net conduction in the y 
direction. This represents the physical condition of a large value of the 
Peclet number. By an analogous reasoning, the radiative heat transfer in 
the X direction can be neglected in comparison to that transferred in the y 
direction. If, in addition, it is assumed that the Eckert number of the 
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Figure 7.1 Physical model and coordinate systems for 
radiating gases between parallel plates. 
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flow is small, then Eq. (7.1) reduces to [8] 


9t 


+ 


3T . 

3T 


u — + V — - 

a 

3X 

9y 

3y2 


pCp 3y 


(7.2) 


where a = (k/pCp) represents the thermal diffusivity of the fluid. 

For a steady fully-developed flow, v = 0, and u is given by the well- 
known parabolic profile as 


u = 6 (5 - ?2); € = y/L (7.3) 

where u^^ represents the mean fluid velocity. Also, for the flow of a per- 
fect gas with uniform wall heat flux, 8T/3x is constant and is given by 

3T/8X = (2aq^)/(u^L/k) (7.4) 

Now, by combining Eqs. (7.2) - (7.4), the energy equation is expressed in 
nondimensional form as 


11 = I2(c - 52 ) = 111 . l_!!i (7.5) 

3t 3?2 9g 

where 

T = at/L2 ; e = (T - Ti)/(q^L/k) 

The expression for 3qj^/9€ in Eq. (7.5) is obtained from either Eq. (3.22) or 
Eq. (3.24). 

By assuming that the initial temperature distribution in the gas is 
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some uniform value Tq = Tj, the initial and boundary conditions for this 
problem can be expressed as 


0(f, 0) = 0 

(7.6a) 

9(0, t) = 9 (1, t) * 0 

(7,6b) 

9^(C = 1/2) = 0 ; 9^(f = 0) s . 0^(5 = 1} 

(7.6c) 


It should be noted that all the boundary conditions given in Eqs. (7.6) are 
not independent and any two convenient conditions can be used to obtain 
solutions. Also, the initial temperature distribution can be any specified 
or calculated value of 9{S, 0) = f(C). 

For flow problems, the quantity of primary interest is the bulk temper- 
ature of the gas, which may be expressed as [41] 

0b “ (Tb-Ti)/(qwL/k) = 6 9{?, T) (? - ?2) df (7.7) 

The heat transfer q^^ is given by the expression, q^ = h^ (Ti-Tj^), where h^ 
is the convective heat transfer coefficient (W/cm^-K). In general, the heat 
transfer results are expressed in terms of the Nusselt number Nu = h^ 

Here, represents the hydraulic diameter, and for the parallel plate geom- 
etry it equals twice the plate separation, i.e., Dpj = 2L. Upon eliminating 

the convective heat transfer coefficient h from the expressions for q and 

c w 

Nu, a relation between the Nusselt number and the bulk temperature is 
obtained as 


Nu = 2 Lq /k(T, -T, ) = -2/e^ 
w 1 b b 


(7.8) 
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The heat transfer results, therefore, can be expressed either in terms of Nu 
or 9^. 

7.1.1 Steady Laminar Flow 

For Steady-State conditions, ae/3t = 0 and Eq. (7.5) becomes 

0" - 12 (5-f2) = (1/q^) dqj^/dc (7.9) 

By integrating Eq. (7.9) once and using the conditions that at c = 1/2, 
qj^(?) and (de/d5) are equal to zero, one obtains 

0 ' - 2(352 - 253 ) + 1 = qR(5)/q^^ (7.10) 

The expression for qR(C) in Eq. (7.10) is obtained from either Eq. (3.21) or 
Eq. (3.23). 

For the present physical problem, e^ =» &2 for 

the case of linearized radiation, a combination of Eqs. (3.21), (5.2a), and 
(7.10) results in 

0 ' - 2 (352 - 253) + 1 

= 1(L/I0 {/« 

- 0(0 C- Uq^ (€'-?)] d5'} (7.11a) 

5 2 

A combination of Eqs. (3.23), (5.2a), and (7.10) gives an alternate form of 
the energy equation for the steady case as 
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0' - 2(3c2 - 2f3) + 1 

" £ - 3 

= (L/k) E H,| r (d0/d5') A.[-u„.U-r)] dS' 

+ (d0/d5') A. [1 u..U'-c)] dC] (7.11b) 

5 1 2 

Note that this equation can be obtained directly by integrating the left- 
hand side of Eq. (7.11a) by parts. Equations (7.11) provide two forms of 
the energy equation for the steady-state conditions. 

For the case of negligible radiation, Eqs. (7.11) reduce to 

0' = 2 (3?2 - 2g3) - 1 (7.12) 

The solution of Eqs. (7.12) is found to be 

0(C) = C (2?2-53.i) (7.13) 

Thus, a combination of Eqs. (7.‘7) and (7,13) gives the result for the bulk 
temperature for the steady case with no radiation as 

- 0,^ - 17/70 (7.14) 

This result is useful in determining the extent of radiative contributions. 

7.1.2 Transient Radiative Interactions 

For the transient case, a combination of Eqs. (3.22), (5.2a), and (7.5) 
gives the energy equation for the linearized radiation as 
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- 8^ - 3N0 - 12 (c-52) 


= -j (L/k) -^0 t-l 


+ f^ e(?', t) AV C-Uq.U'- 5)] d|'} 
C 1 2 


(7.15a) 


where 


N = (PL2/k) K, = (PL2/k) E S.(T) (de, /dT). 

i i*l 1 “i 'i 

Note that this definition of N is slightly different than the definition of 
Nj in Eq. (5.4c). The dimensionless gas property N characterizes the rela- 
tive importance of radiation versus conduction within the gas under opti- 
cally thin conditions [23, 41]. Also, by combining Eqs. (3.24), (5.2a), and 
(7.5) another form of the transient energy equation is obtained as 


- 0^-12 ( 5 - c 2 ) 

= - (L/k) E H,.Uq. rr^(90/35') A1 [- u«.(e-5')] de' 

2 i=l 0 ^ 2 ^ 

- (30/9C) Al [| Uq. (g'-^)] d?') (7.15b) 


Note again that Eq. (7.15b) can be obtained directly by integrating the 
left-hand side of Eq. (7.15a) by parts. Quite often, Eq. (7.15b) is the 
convenient form to use in radiative transfer analyses. 

For the case of negligible radiation, N = 0 and both forms of Eq. 
(7.15) reduce to 
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e?? - ®T " 

By employing the product solution procedure, the solution of Eq. (7,16) can 
be obtained and the result can be expressed in terms of the bulk temperature 
through use of Eq. (7.7). 

The solution of Eq. (7.16) is assumed to be of the form 

e(5,T) = g(c) + h(?,T) (7.17) 

From Eqs. (7.16) and (7.17), there is obtained two separate equations as 

g" = 12(C-?2) (7.18) 

h^^-h^ =0 (7.19) 

The solution of Eq. (7.18) is obtained by direct integration as 

g(C) = U252-C3-1) (7.20) 

This is the same result as given by Eq. (7.13) for the steady case if g(c) 
is replaced by e(?). The solution of Eq. (7.19) is found to be (see Sec. 
7.2) 


n 

h(?,T) = E C„ sin (nirS) exp[-(nir)2T] (7.21a) 

n=l ” 


where 
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C„ = - 2 g(?) sin(n7rS) d5, n = 1,2,... (7.21b) 

"0 

Thus, the complete solution of Eq. (7.16) is given by 


0(C, t) = ? (2^2 -1) 
n 

+ E C„ sin (ac) exp(-a2T) ; a * nir (7.22) 

n»l " 


The expression for is obtained from Eqs. (7.20) and (7.21b) as 


C„ = (4/a5) [12 - 12a2 + a*») cos(a) - 24] , n = 1,2,.... (7.23) 

where a is defined in Eq. (7.22). By combining Eqs. (7.7) and (7.22), the 
expression for the bulk temperature is obtained as 

00 

= -17/70 + 6 E C [(1/a) + (4/a3)] exp (-a2T) (7.24) 

° n=l " 


where is given by Eq. (7.23). 

7.2 Op t i cal 1 y Thin Li m i t 

In the optically thin limit, the steady-state energy equations, Eqs. 
(7.11a) and (7.11b), reduce to 

6' - 2(3c 2-2?3) + 1 = i N[/^ e(C') dC' - 0(5') dC'] (7.25a) 

2 0 0 


e' - 2 ( 3 ^ 2 - 253 )+ 1 = i N [/^ (?-5‘) (de/d5')d5' 

2 0 
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+ (g '-c) (d0/d5')df (7.25b) 

The differentiation of Eqs. (7.25a) and (7.25b) yields the same energy 
equation for the optically thin limit as 

9" - 3N0 = 12 (c-?2) (7.26) 

The solution of Eq. (7.26) satisfying the boundary conditions 9(0,t) = 0 and 
0(1,t) = 0 is found to be 

0(C) = (16/3N2) [sinh(-V^2)/sinh (V^)l coshi^lN'(c-l/2)] 

+ (4/N)(?2-C + 2/3N) (7.27a) 

Alternately, the solution of Eq. (7.26) is written as 

0(C) = Cl exp (V^C) + C 2 exp (- ?) 

+ (l/m2 ) (24-12 me + 12 m 52 ) ; m = 3N 

The constants Ci and C 2 are obtained by using the boundary conditions ®(0) = 
0 and 0'(l/2) = 0, and the solution for 0(c) is found to be 

0(C) = (l/m2) {-24/(l+e‘V^] (e'V^eV^ ^ + e‘V^ ^) 

+ 24 - 12 m C + 12 m C^} (7.27b) 

Equations (7.27a) and (7.27b) should produce identical results. The 
expression for the bulk temperature, in this case, is obtained by combining 
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Eqs. (7.7) and (7.27b) as 

«b = ^ Ca-e’V^/d+e'V”)] - iE + li - ii (7.28a) 

ra7/2 m3 m2 5m 

or 

- 1/2 

= [1/(3N)3] {576(3N) (NEXP) - 21.6N2 + 72N - 288} (7.28b) 

where 

1/2 1/2 
NEXP = {1-exp [-(3N) ]}/{l + exp[-(3N) ]} 

In the optically thin limit, the transient energy equations, Eqs. 
(7.15a) and (7.15b), reduce to 

9^^-0 - 3Ne = 12 (?-c2) (7.29a) 

T 

e - 0 - 12(5 - 52 ) 

55 T 

= I n [f^ (30/05 ')d5'- / (36/35') d5'] (7.29b) 

2 0 5 

Note that Eq. (7.29b) is identical to Eq. (7.29a). 

The solution of Eq. (7.29) is assimied to be of the form 

0(5,t) = g(5) + h (5,t) (7.30) 

Thus, Eq. (7.29) can be written as 

hgg - \ - 3Nh = - g^g + 3Ng + 12(5 - C^) (7.31) 

Consequently, 
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g" - 3Ng = 12(? - (7.32) 

and 

^5? - \ - 3Nh = 0 (7.33) 

N 

The conditions for Eqs. (7.32) and (7,33) are obtained from Eq. (7.6) as 

0(O,t) = h(0,T) + g(0) * 0; h(0,T) = 0, g(0) = 0 (7.34a) 

9(1,t) ^ h(l,T) + g(l) = 0; h(l,t) = 0, g(l) = 0 (7.34b) 

0(5,0) = h(5,0) + g(5) - 0; h(5 ,0) = - g(5) (7.34c) 

The solution of Eq. (7.32) satisfying the boundary conditions given by 
Eqs. (7.34a) and (7.34b) is identical to the solution of Eq. (7.26) as given 
by Eq. (7.27) if 0(5) is replaced by g(5), i.e. , 

g(5) = (16/3N2) [sinh (- ^^iN/2)/ sinh(V3N)] cosh [VsN (5-1/2)] 

+ (4/N) (52 - 5 + 2/3N) (7.35) 

The solution of Eq. (7.33) is obtained by using the product solution 
procedure and implying the conditions h(0,T) = h(l, t) = 0 and h(5,0) = 
-g(5). For the product solution, it is assumed that 

h(5,T) * F(5) G(t) (7.36) 

By using Eq. (7.36), Eq. (7.33) is separated into two ordinary differential 
equations which are expressed along with appropriate conditions as 


F" + X2 F = 0 ; F(0) = 0, F(l) = 0 


(7.37) 
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6 + (3N + x2) G = 0 ; h(5,0) = 6(0) = g(c) (7.38) 


The solution of Eq. (7.37) is given by 


F^(5) ~ sin (nir^)* n ~ 


and the solution of Eq. (7.38) is found to be 


6^(t) » exp {- [3N + (mr)2] x} 


Thus, the complete solution of Eq. (7.33) is 


where 


h(5,T) = l sin (nirC) exp {-[3N + (nTr)2]T} 
n=l " 


= -2 g(c) sin (nir?) d?, n = 1,2,... 
" 0 


(7.39) 


(7.40) 


Now, the solution of Eq. (7.29), as expressed by Eq. (7.30), is written as 


9(C>t) = (16/3N2) [sinh (-^Ih/Z)/ sinh(^ 3N)] cosh (5-1/2)] 

+ (4/N)(52 - 5 + 2/3N) 

CO 

+ s sin (nir5) exp{-[3N + (nTr)2]T} (7.41) 

n=l " 

From Eqs. (7.35) and (7.40), it follows that 
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“ 0 , for n even (7.42a) 

= 32C3N+(mr)2]/[3N2(mr)3] +2(nTr)/{3N2[3N+(nir)2]), for n odd (7.42b) 


By combining Eqs. (7.7) and (7.41), the expression for the bulk 
temperature is obtained as 

= 6 { (16/3N2)[sinh(-Vii/2)/Sinh(/ii)][(l/^^ cosh( V5/2) 

- (4 +V^) (3N)-3/2 sinh(Vw2)] + (4/N)[-l/30 + 1/(9N)] 

00 

+ s C [l/(nir) + 4/(nir)3] exp[-(3N + n2ff2)T]} (7.43) 

n=l " 

where is given by Eq. (7.42). 

7.3 Large Path Length Limit 

In the large path length limit, the steady-state energy equations, Eqs. 
(7.11a) and (7.11b), reduce to 

0' - 2 ( 3^2 - 2c 3) + 1 = M /' 0(5') d5'/(? - 5') (7.44a) 

0 

0' - 2(3c2 - 253) + 1 

" p 1 

= (L/k) Z H,.{/Md0/d5') tn[i. u..(5 - 5')]d?' 
i=l ^^0 2 

+ /^(de/dc') in [I uq^U' - 5 )] d 5 ' (7.44b) 

where 

n 

M = Hi L/k = (L/k) AQ.(de^./dT)^^ (7.44c) 

Through integration by parts, it can be shown that Eq. (7.44b) reduces to 
Eq. (7.44a). The parameter M in Eq. (7.44c) is defined differently than 
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In Eq. (5.4e) . TTie nondimensional paraneter M constitutes the radiation- 
conduction interaction parameter for the large path length limit [23, 41]. 
Equation (7.44a) does not appear to possess a closed form solution; a 
numerican solution, however, can be obtained easily. 

In the large path length limit, the transient energy equations, Eqs. 
(7.15a) and (7.15b), reduce to 

“ 12 (? - ? 2 ) 

= - (HiL/k) [ /S(g',T) dc'/(? - t')2 + /J 6 ( 5 ',t) dc*/( 5 ' - 5 ) 2 ] (7.45a) 

- 0 ^ - 12 (c - ? 2 ) = (HiL/k) (36/3? ')d€ 7(5 - 5') (7.45b) 

Since |(?-5 ')^ | = |(5*-C)^| , Eq. (7.45a) can be written as 

0 ^^ - 0 - 3N0 - 12(? - ? 2 ) 

55 T ^ 

= - (HiL/k) / e(?',T) d5V/(5 - 5')2 (7.45c) 

0 

Through integration by parts, Eq. (7.45c) can be expressed as 
0 __ - 0 ^ - 3N0 - 12 (? - ?2) = (HiL/k) (30/3?') d?'/(? - ?') (7.45d) 

I 

Equations (7.45a) - (7.45d) represent different forms of the governing equa- 
tions in the large path length limit. With the exception of the term (-3N0) 
on the left-hand side, Eq. (7.45d) is identical to Eq. (7.45b). Since N 
represents the radiation-conduction interaction parameter only in the 
optically thin limit [23], it should not appear in the governing equation 
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for the large path length limit. Thus, Eq. (7.45b) is the correct equation 
to use for solution in the large path length limit; the solution of this 
equation is obtained by numerical techniques. 

7.4 Method of Solution 

The solution procedures for both steady and unsteady cases are pre- 
sented in this section. In principle, the same numerical procedure applies 
to both the general and large path length limit cases. 

7.4.1 Steady- State Solutions 

The general solution of Eq. (7.11a) or Eq. (7.11b) is obtained numer- 
ically by employing the method of variation of parameters. For this, a 
polynomial form for 6(5) is assumed in powers of c as 

0(5) = S a s'" (7.46) 

m=0 

By considering a five term series solution (a quartic solution in 5) and 
satisfying the boundary conditions e(0)= 9 '(1/2) = 0 and 0'(O) = -e'(l), 
one obtains 


0(5) = ai(5 - 253 + CM+ 32(52 - 25^ + (7.47) 

Thus, 

e'(5) = 3 i(l - 652 + 4|3) + 32 (25 - 652 + 4g3) (7.48) 

A substitution of Eq. (7.48) in Eq. (7.11a) results in 


ai (l - 652 + 453) + 32 (25 - 552 + 45^) - 2(352 _ 2 g 3 ) + 1 
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= -(L/k) Hi.Uo^{/^e{r) AJ[iuQ.U -5-)]de- 

-/^(c')7r; [iuQ.U' -?)]d5'} ( 7 . 49 ) 

5 1 2 


where expressions for 0(?') are obtained from Eq. (7.47). 

The two unknown constants a^ and ag in Eq. (7.49) are evaluated by 
satisfying the integral equation at two convenient locations (?=0 and 5 = 
1/4 in the present case). The entire procedure for obtaining a^and ^2 is 
described in Appendix E from which it follows that 

ai = (1102 - 16 04 ) /DEN 
a 2 ” (16 03 “ lloj)/DEN 

where 

DEN = 16 (ojo^ “ 0203 ) 

and coefficients oi through 04 are defined in Appendix E. 

Now, with known values of a^ and a2» Eq. (7.47) provides the general 
solution for e(?). The expression for the bulk temperature is obtained by 
combining Eqs. (7.7) and (7.47) as 

ei 3 = (1/70) (17ai + 3 a 2 ) (7.51) 


(7.50a) 

(7.50b) 

(7.50c) 


Note that for the case of no radiative interaction 02* 03, and 04 are equal 
to zero and oj = 1. Thus, aj = 0 and a = -1, and Eq. (7.51) gives the 
result of Eq. (7.14). 

The governing equation for the large path length limit is Eq. (7.44a). 
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For this equation also the solution is given by Eq. (7.47) but the values of 
a's are completely different in this case. There are two approaches to 
obtain solutions in the large path length limit. One approach is to make 
use of Eq. (7.44a) and go through the entire numerical procedure described 
in Appendix E for the general solution. Another approach is to work with 
the general solution but evaluate all and S^. integrals of Appendix E in 
the large path limit. In the large path length limit, the integrals can be 
evaluated in closed forms. Both procedures are described briefly in Appen- 
dix E. In order to distinguish the large path length limit solution from 
the general solution, constants a^ and aj are replaced by bj and bj, and 
coefficients a\ through oti, are replaced by through $4 . The solution for 
the large path length limit, therefore, is given by 


9 (C ) = bi (5 -253+54 ) + b2 (C2 -253 4^*^ ) 

where 

bi = (11b 2 - I 6 B 4 ) /BOTTOM 
b 2 = (1663 - lie 1 ) /BOTTOM 
BOTTOM- 16(BiB4 - 92 ^ 3 ) 


(7.52a) 


(7.52b) 

(7.52c) 

(7.52d) 


and coefficients gj through 04 are defined in Appendix E. For this case, 
the bulk temperature is given by 


= (1/70) (17bi + Bbj) 


(7.53) 


Note that in this case the value of coeffients 61 through 64 are obtained in 
closed form. 
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7.4.2 Transient Solutions 

The governing energy equations for the transient case are Eqs. (7.15a) 
and (7.15b). As in Sec. 6.3, general solutions of these equations are 
obtained numerically by employing the method of variation of parameters. 

For the present problem, a polynomial form for e(c,t) is assumed as 

0(?,t)-E an,(T) l"* (7.54) 

m=0 


For a quadratic temperature distribution in 5 (with time dependent coeffi- 
cients), Eq. (7.54) is written as 

e(?,t) = ao(T) + ai(T) 5 + a 2 (T) (7.55a) 

By using the boundary conditions 0(O,t) = 0 and 9^(?-l/2) = 0, this reduces 
to 

9(5.t) = g(T) (5-?2) (7.55b) 

where g(t) represents the time dependent coefficient. Consequently, 

0g(?.T) = g(T) (1-25); = -29(t); 0^(5,t) = (5-52)g'(T) (7.56) 

Also, a combination of Eq. (7.6a) and (7.55b) yields the initial condition 

e(C,0) = g(0) = 0 (7.57) 


Note that essential boundary conditions are used already in obtaining the 
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solution represented by Eq, (7,55b). 

By employing Eqs. (7.55b) and (7.56), Eqs. (7.15a) and (7.15b) are 
transformed in alternate forms which are expressed in a compact form as 


9'(t) + 


0i(5); 

J2(5). 


9(t) + 12 = 0 


(7.58) 


where Ji(c) and O 2 U) af*e defined in Appendix F. Hie function Ji(5) is used 
for solution of Eq. (7.15a) and 82 ( 5 ) is used for solution of Eq. (7.15b). 
The solution of Eq. (7.58) satisfying the initial conditions of Eq. (7.57) 
is given by 

g(x) » J£-{exp [-J(5)t]-1) (7.59) 

0 («) 

The temperature distribution given by Eq. (7.56b) can be expressed now as 


e(?,r) » Ji-{exp (7.60) 

0(5) 

The expression for the bulk temperature is obtained through use of Eq. (7.7) 
as 


0 = 72 /^[(C-c2)Vj(5)]{exp[-J(5)T]-l} (7.61) 

o 0 

Note that in Eqs. (7.59)-(7.61), J(?) becomes Ji(5) for solution of Eq. 
(15a) and 82 ( 5 ) solution of Eq. (15b). 

For a quartic solution in C, Eq. .(7.54) gives the result identical to 
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Eq. (7.47) which for the transient case is expressed as 


9{5.t) = g(T)(C - 253 + ?*») + h(T)(c2 - 253 + s'*) (7.62) 

Thus, 

9g (S,t) = (1 - 652 + 453)g(x) + 2(S - 3S2 + 2s3)h(T) (7.63a) 

955(5.0 = 12( - £ + £2) g(x) + 2(1 - 6C + 6£2) h(x) (7.63b) 

0^.(5 .f) = (£ - 2£3 + ?*») g'(T) + (?2 - 2£3 + £‘*)h'(T) (7.63c) 


By substituting Eqs. (7.62) and (7.63) into Eq. (7.15a), one obtains 
X g'(x) + J3(£) g(x) + yh' + A(£) h(T) = - z (7.64) 

where 

X = (c - 253 +?*» ); y * (€2 - 2g3 + ?4 ). z = i2(C - 52 ) 

and functions 33(5) and <!,(£) are defined in Appendix F. Equation (7.64) 
constitutes one equation in two unknowns, namely g(x) and h(x). However, 
since the equation is linear in x, the principle of superposition can be 
used to split the solution into two solutions as 

X g'(x) + J3(£) g(t) = -z/2 (7.65) 

yhV(x) + h(x) = -z/2 (7.66) 

The initial condition for this case can be written as 

9(£.0) = g(0) (£ - 2£ + £3) + h(0)(£2 - 2£3 + 5*^) = 0 (7.67a) 

Consequently, 


g(0) = 0; h(0) = 0 


(7.67b) 
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The solution of Eqs. (7.65) and (7.66) satisfying the appropriate 
initial condition of Eq. (7.67b) is given respectively as 

g(t) = [z(5)/2J3(5)] (exp[-03(5)T/x(5)]-l) (7-68) 

h(r) • [z(C)/2J,(5)] (exp[-J^(|),/y(5)]-l) (7.69) 

By substituting Eqs. (7.68) and (7.69) intp Eq. (7.62), the expression for 
the temperature distribution 1s obtained as 


9(5,r) ' [6U-?7)(6-2sZ+{S)/J3({)](exp[-03(?)r/x(? )]-!} 

+ [6(5-57)(?7-2E3a«'.)/J^(|)](exp[-J^(5)T/y(5)]-l} (7.70) 

The bulk temperature in this case is given by 

9^ . 36 [(5 - tZ)(| - 2?3 + 5'*)/03(5)](exp[-J3(5)T/x(5)] -lid? 

+ 36 /[(6 - S2)(|7 - 2e3 + 5M/04(5)]{exp[-Js(5)x/y(l)]-l}d5 (7.71) 

where x and y are defined in Eq. (7.64). 

By substituting Eqs. (7.62) and (7.63) into Eq. (7.15b), there is 
obtained 


xg' + JsCc ) g(T) + yh' + J6(5) = -z (7.72) 

where again x,y,z are defined in Eq. (7.64) and functions JsC?) and Je(^) 
are defined in Appendix F. The solution procedure for this equation is 
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identical to that for Eq. (7.64) and the results for temperature distri- 
bution and bulk temperature are given respectively by Eqs. (7.70) and (7.71) 
with Jg replaced by Jg and by Jg. 

In the large path length limit, the two applicable governing equations 
are Eqs. (7.45b) and (7.45d). The solutions of these equations can be 
obtained from the general solutions by evaluating the integrals in J 
function in the large path length limit. 

Alternately, for a quadratic temperature distribution, Eqs. (7.45d) and 
(7.45b) are transformed respectively to 


g'(T) + 


J7(5) 

08 ( 5 )’ 


g(x) + 12 = 0 


(7.73) 


where 07 ( 5 ) and <^(5) are defined in Appendix F. The solution of Eq. 
(7.73) is given by Eq. (7.59) and expressions for 0(c,t) and can be 

obtained from Eqs^ (7.60) and (7.61) respectively. Of course, proper care 
should be taken to use the correct relation for 0 functions for different 
equations. The large path length limit solutions for a quartic temperature 
distribution can be obtained in a similar manner. 



8. PLANS FOR SPECIFie RESULTS 


Some specific results have been obtained and these are being analyzed. 
The present plans are to obtain extensive results for the following cases with 
varying physical and flow conditions: 


A. Physical Geometries 

1. Parallel Plates: One-Dimensional Radiation 

2. Parallel Plates: Two-Dimensonnal Radiation 

3. Diffusing Channel Flow: One- and Two-Dimensional Radiation 

4. Channel Flow; Top Plate Flat, Bottom Plate with a S-IS" ramp 

5. Scramjet Inlet Configurations 

B. Radiative Interaction Cases 

1. Transient Radiative Transfer in Homogeneous Gases 

2. Transient Energy Transfer By Radiation and Conduction in Homogeneous 
Gaseous Systems 

3. Transient Energy Transfer By Radiation, Conduction, and Convection in 
Homogeneous Gaseous Systems 

4. Applications to Flow of Homogeneous Gaseous Mixture 

5. Applications to Flow of Chemically Reacting Gaseous Mixtures 

6. Applications to the Scramjet Inlet Configurations. 

C. Boundary Conditions 

1. Isothermal Black Boundaries 

2. Isothermal Gray Boundaries 

3. Nonisothermal Boundaries 

4. Boundaries with Uniform Heat Flux 

5. Actual Scramjet Inlet Configurations 

D. Flow Conditions 

1. Incompressible - Various Cases 

2. Compressible - M^ = 0.5, 1, 2, 3, 4, and 5 

3. Temperature Range - 300, 500, 1000, 2000, and 5000 K 

4. Pressure Range - 0.1, 1, 2, 5, and 10 atm. 

5. Realistic Conditions for Scramjet Inlet Flows 


E . Pa r ti c i pa ti ng Me d i urn 

1. CO - One Fundamental Band 

2. COo - Three Important Bands 

3. HpO - Five Important Bands 

4. OH - One Fundamental Band 

5. CO 2 + HpO (Different Concentrations) 

6. OH + H 2 O (Different Concentrations) 

7. OH + H 2 + O 2 + H 2 O (Different Coneentrations) 

8. A Realistic Combustion Model for the Hydrogen Air System 



Specific Results 


1. Nongray Solutions Based on Band-Model Correlations 

2. Optically Thin Solutions 

3. Large Path Length Limit Solutions 

4. Gray Solutions 

5. Nongray Solutions for Scram jet Inlet Flows 



9. CONCLUDING REMARKS 


A brief review is presented on various band models and band model 
correlations that are useful in nongray radiative transfer analyses. 
Different formulations for one-dimensional radiative flux are provided. These 
are used to develop the basic governing equations for transient energy 
transfer in gaseous systems. Limiting forms of these equations are obtained 
in the optically thin and large path length limits. Numerical procedures are 
described to solve the governing equations for different physical and flow 
conditions. The plans for obtaining extensive results for different cases are 
provided. The formulation and numerical procedure presented in this study can 
be extended easily to multi-dimensional analyses. In the near future, the 
influence of radiative interactions will be investigated for the realistic 
flow of hydrogen-air mixture in the scramjet inlet configuration. 
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APPENDIX A 

EXPONENTIAL INTEGRALS AND EXPRESSIONS FOR RADIATIVE FLUX 
Some important relations for the exponential integrals are given 
Appendix B of Ref. 8; and it is noted that 


/ E„(t) dt = - E„ ^ j (t); E„(0) • l/(n-l) . 

Now, consider the first integral in Eq. (3.1) as 

UD = fl’' Ej (tj.t)'’‘ 

By defining x = t - t, dt * - dx, and Eq. (A. 2) becomes 

A 

1(1)=/° E2(x)(-dr) = ;’' E2(x)dx 

= - [E3(x)]/ = E3(0) - E3 (Tj) . 

Thus, 

E3(t>) = £3(0) - 1(1) = i-/? E3(xj-t)dt 
The second integral in Eq. (3.2) is written as 


in 


(A.l) 


(A.2) 


(A.3) 


1(2) E2(t-T2)dt . 


(A.4) 
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8y defining x = t - Eq. (A. 4) is expressed as 
U2) = E 2 (x)dx = [E 3 (x)I 

' - rE3<"0X-‘x> - ^3(0)] 
or 

EgCxox-Tx) = £ 3 ( 0 ) - 1(2) =i - E2(t-T^)dt . (A.5) 

2 A 

A substitution of Eqs. (A. 3) and (A.5) into Eq. (3.1) results in (for 
the case = ©ix and B 2 X = e 2 x) 


’^RA ®1a - ®2x 

+ 2 {/q^ (t) - eiA] E2(Tj^-t) dt} 

- 2 [%x(t) - e 2 X]Ej(t-T,)dt} (A.6) 

X 

This equation when expressed in terms of the wave nimiber oj is exactly the 
same as Eq. (3.4). Following a similar procedure, Eq. (3.5) is obtained 
from Eq. (3.2). 
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APPENDIX B 

ALTERNATE FORMS OF RADIATIVE FLUX EQUATiONS 
In radiative formulations, it is desirable to express the relations for 

— —I 

and div q^ in terms of A and A , and avoid the use of A . This 
r r 

, I II 

is accomplished by expressing the integrals containing A and 7T in 
alternate forms through the procedure of integration by parts. This is 
performed by using the relation 


/, m dn = (mn)^ - n dm 

a u a 


(B.l) 


Consider now the first integral in Eq. (3.21) and express as 

BI(l) = fl (?•) a'[| Ug^. (5-?V)] d?' (B.2) 

For integration by parts, let 


Then, 


m = Fiu; dn 


a' [i u (c-5V)l dC . 
2 ° 


dm = (dFiu/d?') dC - (de/d^') d?' . 
In order to get n, let 

u =1 uq du/dC = - 3 Uq/2 


and 

dn = [dA(u)/du] d5' = [dA{u)/dC' x dc'/du] de' 

= - (2/3Uq) [d A(u)/d5'] dV = (- 2/3uQ)dA(u) . 
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Thus, 

n = -(2/3 Uq) A(u) = -(2/3Uq)A [| Uq(?-C')] • 

Consequently Eq. (B.2) can be written as 

BI(l) = Fia,U)r- — A(0)] - Fio,(0) J (1 U. 5 )] 

3u„ ^ 3u„ 2 ° 

0 0 

+ — /„ [de (5')/d5'] 1[ iu d-5')]d{' (B.3) 

3Uq ° ^ ° 

Note that by the definitions given in Eqs. (2.2) and (2.8), A(0) = 0. In 
the present case, only the definition given in Eq. (2.8) is acceptable for 
A(0) = 0. Also by definition, FjojCO) = 0. Thus, Eq. (B.3) reduces to 

Bid) = m%) !l [de^(5')/de1 A [| (8.4) 

llie second integral in Eq. (3.21) is written as 

BI(2) = (2/3Uq) fl IdejV) A [I u^ (B.5) 

m = F 2 U; dn = A'[j UQ(c'-5)]dc' 


Let, 
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Then, 


dm = ( dp2^/d5 ' )d? ' = { de^/d? ' ) d5 ' 
Now, to get n, let 

u = ^ Uo(C ' , du/d? ' = 3 Uq/2 


and 

dn = [dA(u)/du]d? ' = [dA(u)/d?' x d?'/du]d?' 
= (2/3uQ)[dA(u)/d?'i d?' = (2/3 Uq) dA(u) 


Thus, 

n = (2/3 Uq)A(u) =(2/3Uq)A[1 Uq(?'-?)] 

Consequently, Eg. (B.5) is expressed as 

BI(2) = F2„(1) {(2/3Uj,)A [iu^(l-5)])-{F2,.U) [(2/3u^) A(0)]} 

- (2/3Uj,) l\ [dej5')/ds'] A [1 11^(5'-?)] d{' (B.6) 

Since, by definition, ^ reduces to 
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BI(2) = - (2/3u.) [de (?')/dC'] A [1 u (§'-?)] dC (B.7) 

U ^ CD 2 ^ 

By use of Eqs. (B.4) and (B.7), the integrals in Eq. (3.21) can be 
expressed in alternate forms and this results in Eq. (3.23). 

Consider now the integrals in Eq. (3.22); the first integral is written 
as 

BI(3) = fl Fi^(C‘) a' [1 Uq (C-s')] d?' . (B.8) 

For integration by parts, let 

m = Fjo) (s'); dn = A [j u^ (?-?')] d? ' 

Then 

dm = (dFi^/d^') d^' = (de^/df ) d^' 

As before, to obtain n, let 

u = |- Uq(s-S'), du/ds' = - 3 Uq/2 


and, therefore. 
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dn = ^ f d7r(u)/du] d^' 
du 

= {_! ^’fdA(u)/du] dg' 
dr du 

= df(-2/3 Ujj)fdA(u)/du]} 

Thus, 

n = (-2/3 Uq) a'(u) = (-2/3 u^) A [i u^ (€-r )] 

and Eq. (B.8) can be written as 

BI(3) • {Fi» (-2/3 uj a' (5-5')]lo 

-/o((-2/3 %) a'[|u^(C- 5')] )[dF,»(?')/dS'] dV- 

I 

Since FiuCO) = 0 and A (0) = 1, Eq.(B.9) reduces to 

61(3) = (2/3u^) (-F,„(5) * fl [de^(c')/dc'] A'[l (C-5')]d?') 

Similarly, for the second integral in Eq. (3.22), one can find 
BI(4) =/‘ Fj.. (5') a'[|u^ (5'-{)] dE' 


( 8 . 9 ) 


(B.IO) 


= (-2/3u^){F2(e) + [de_^ (5’)/d{'] A'[l (5‘-«)]d5'} (6.11) 
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A substitution of Eqs. (B.IO) and (B.ll) in Eq. (3.22) results 
(3.24). Also, a differentiation of Eq. (3.23) with respect to C, 
the Leibnitz formula, gives 


dq,.U) n .de^i (5) 


dc i=l 




d? 


7^(0) - 0 


3u 


. C de .(5') , ^ 




+ 0 . ^ilL!ix.(o) 

d5 ' 


3u 


01 


de.(C') 


(Ul 


dC 




Since A(0) = 0, Eq. (B.12) reduces to Eq. (3.24). 


in Eq, 
by using 


(B.12) 
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APPENDIX C 

DEFINITION AND EVALUATION OF INTEGRAL FUNCTIONS 
For the convenience and use in the computationaT procedure, the follow- 
Ing definitions are employed in expressing the relations for the integral 
functions: 


b . = 3u 72 
1 01 

c, = 1/b, = 2/3u^, 

C(5) = 1/(5 -SM 
r = 1 /(ot|) 

Various integrals are defined and simplified as follows: 

Gi(5) = C(?) (3Ni(5-52) + 

-I- i . u2 . {/^ (c • -? ) 7?. " [b. {? -C ' )] d5 ' 

(5 • -€ '2 ) A^. " [ b^ (C • -C ) ] dC ' } ) ( C . 1 a) 

= C(C) {3Ni (5-C2) + 

+ 1 Uq. {/g^ [c-C|U-c2+2cc^u-(c^u)2] J (u)du 

bj(l-5) 2 

+ / ^ [?-»-c^.u- 52 . 2 cc^.u-(c^ u) ]‘A''(u)du}) 

0 


(C.lb) 
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W = C(5) (3Ni(i- 52 ) 

■? ii "oi i/o Ai"[b^(5-CV)]d?' 

+ (5'^-l) a! [b.(5'-?)]dC'}) 

= CU) (3Ni (1 - 52 ) 

2 

- I I "li U0I <V' ^ 

bi(l-S) 

+ !q [ 5 ^ + 25 e. u + (c^u)2 -1] I. (u) du}) 

63(0 = C(5) (i I {/; (1-25 •) A^ [b^ (5-5 ' jldg ' 

2 r=Z ‘ 

- (1-25') [b^ (5'-5)]d5'}) 

= C( 5 ) (_?_ ^ [ 1-25 + 2 c^u] J.(u)du 

b.(l-5) 

- Jo [1-25 - c^u] A^.(u)du}) 

" p _' 

a.(5) = C(5) (-3 2 {/q 5' A^ [b.(C4-)]dr' 

i *1 

- /| c a; [b, (5'-5)]br)) 


(c. 2 a) 


C.2b) 


(C.3a) 


(G.3b) 


(C.4a) 
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» C({) C- 2 I M [ 5 -c u] 7 ?! (u) du 

1=1 ' ' 

b^(l- 5 ) 

+ [5 + c^u] A!(u)du}) 

65(5) = { 3 r/ 8 ) Uq,. Hj. ill a! Ib^{ 5 -r)]d€' 

n b.? 

= (r/ 4 ) [5-c^u -?2 + 2 ?c.u-(c^u) ] Aj(u) du 

- /q [«+c^u-c2 - z^c.u -(c.u) ] X (u) du} 

Ge (0 ? (3 r/8) Uq^. Hj. Ul A! [b^. ( 1 - 5 ' )]d 5 ' 

+ /g ( 1 - 5 '^) j! [b^ (5'-5)] d?'} 

n b 5 

= (r/ 4 ) [52-25 c^u +(g^u) 2] I. (u) du 

b 4(1-5) 

+ /q [1-52-25 c.u - c>u) 2 ] A.(u) du} 


(C 4 .b) 


(C. 5 a) 


(C. 5 b) 


(C. 6 a) 


(C. 6 b) 
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G7(c) = (r/4) {fl d-^') [b^(C4')] 6V 

-f^ (1-20 fb^(f-C)] d5'} 

= (r/6) (Hi^/Uq^) (1-2c + 2e.u) 7r.(u) du 

b.Cl-c) 

- (1-2s-2c^u) A^. (u) du} 

GsCn = (r/4) ill (20 A^ [b,.(?-5 ') d5 • 

- /^ (20 7?^ [b. (5 '-c)] dc'l 
C 

n b .5 

= (r/6) r ^(Hii/uQ^){/jj (2c - 2 g^u) (u) du 

bdl-€) 

- /q (2c + 2c. u) J^iu) du} 

G, = (3T/3) u„i Hj, {/‘ (5'.5'2) A, (b,5 ' ) dC) 

n b. 

= (r/4) [c.u - (C|U)2] A^ (u) du} 


(C.7a) 

(C.7b) 

(C.8a) 

(C.8b) 

(C.9a) 


(e.9b) 
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Gio = (3r/8) Uoi {/q (1-5 ') Ai (b^5' ) dS'} 

n 1^ 

= (r/4) Hj. {/g'' (1-c.u) J- (u) du} 

n 1 

Gil = (r/4) {/q (1-25') (b^S') d?'} 

n li) 

= (r/6) (H^. (1-2 c^u) 7T. (u) du} 

1 _ 

Gi 2 » (r/4) .^(Hii/Uq^) {/q (25 •) A^ (b^ 5') d5'} 

n ^ 

= (r/6) .I^(Hi^/Uo^.) {/q’ (2 c> u) A. (u) du) 



(C.lOa) 

(C.lOb) 

(C.lla) 

(C.llb) 

(C.12a) 

(C.12b) 

(C.13) 


(improper) 


(C.14) 
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8,5 (5) = Ml C(5) /' d-X') -il- (C.lSa) 

( 5-0 

= Ml C(5) {2 + (25-1) [(5-l)/5]} (C.15b) 

G„ ({) ■ -2M, C({) s ' . '* ■ 1' — (C.16a) 

(5-5 ') 

= 2 Ml C(5) {1+ 5 An [(5-l)/5]} (C.16b) 

G,, ({) • (rH,/4) ({'-c'2) (C.17a) 

= (rH,/4) { - 1 + { - ({ - 52 ) tn [(e-D/ 5 ]) (C.17b) 

Gi 8 (?) = (rHi/4) rr’ e '2 —ill- /’ — ] (Improper) (C.13) 

(1-5') ( (?-«') 

n 

Gi 9 (5) = (r/4) E H,. fr (1-25') An [b4(5-5 ' )]d5' 
i=l 0 

+ /J (1-2?') An [b. (?'-?)] d?'} (C.19a) 

* (rHi/4) {- i + ? + (?2-?) An [(1-?)/?]} (C.19b) 


G20 (5) = (r/ 4 ) Hii Ul 25' An [bi(5-5')] d5' 



78 


+ 2C tn [b. (c'-e)] dc'} 

- (r/4) (H ,{52 tn U/(l-c)3 - 5 - 

+ X H,, Jln [b,a-e)3) 
i=l ^ ' 

G 21 = (rHi/4) /o U'-5’2) dc7C = (rHi/8) 

G 22 = (rHi/4) /q (1-5 '2) d5'/f 

= - (rHj/4) [i + in (0)] (improper) 

n 

G 23 = (r/ 4 ) I H,. L (1-25 •) Zn (b. ?') d5' 

i_l 

= - (r Hi/8) = - Hi/(8aTi^) 

1 

624 = (r/4) /q 25' iln (b^ 5') d5' 

= (r/4) [- (Hi/2) + I An (b/)] 

i=l 


(C.20a) 


(C.20b) 

(C.21) 

(C.22a) 

(C.22b) 

(C.23a) 

(C.23b) 

(C.24a) 


(C.24b) 
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APPENDIX D 


1. Data: 


INFORMATION FOR NUMERICAL PROCEDURE 


T = 300, 500, 1000, 2000, 5000, etc. ~ K 

P = 0.1, 1, 2, 5, 10, etc. ~ atm 

L = 1, 2, 5, 10, 20, etc. ~ cm 

5 X = 0.0‘-»- 1.0; X = U(I) 


2. Thermal conductivity of the gas: 


Kf = (T, P) 

3. The Planck function and its derivative: 


(D.l) 


'bto 


(T) = 


Cl 

exp (Cs/T) - 1 


= PF 


(D.2) 


'^®btu _ ^2 exp (C 2 /T) 

dT T2 [exp (C 2 /T)-iI^ 


(0.3) 


where 

Cj = (2irhc2)(ia^ 

C 2 = (hc/k) M 
Cj C 2 = (2tt h 2 c3/k)co'* 

and u) = 0 )^ = 0 )^ = wave number center of the ith band, 
c 0 
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By defining 

TEXP = expCCj/T) 

Eq. (0.3) is expressed as 

PFD = (CiC2)(TEXP)/C(T2)(TEXP - 1)2] (D.4) 

All values in Eqs. (0.2) - (0.4) should be evaluated for each band. 

4. Information or relation for for each band, A^ = f(T). 

5. Information on C 2 for each band, C 2 = f(T). 

6. Information on u^ for each band, u^ = 0^2 PL (nondimensional). 

7. Information on B2(T) for each band (nondimensional). 

8. Equivalent or effective pressure relation for each band. 

Pgi “ ^^i ■ (nondimensional) . (0.5) 

In Eq. (0.5), P^ - 1 atm and, therefore, P^ and P must be in the units 
of atm. Note that P. is the partial pressure of the ith species in a 
gases mixture and P is the total pressure; b. (the self broadening 
coefficent) and n are different for different bands. For a single 
component systan, Eq. (0.5) is usually expressed as 

P„ = (b P/P-)'^ (nondimensional) (D.6) 

C U 


9. Line structure parameter for each band: 
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e = 82 Pg = BETA (D.7) 

10. Correlation for each band (for Tien and Lowder's correlation): 

f(8) = 2.94 [1 - exp(-2.60e)] s F (D.8) 

11. Band absorptance correlation for each band (Tien and Lowder's 
correlation) : 

J (u,8) = in [uF + 1] = AU (0.9) 

u + 2F 

12. The derivative of the band absorptance correlation for each band (Tien 
and Lowder's correlation) 

A '(u,6) = [F (u2 + 4 u F + 4F)]/DEN = AUD (D.IO) 

where 

DEN = (F (u2 + 2u + 2) + u](u + 2F) 

These basic relations are used in the governing equations of Section 6 


to obtain nunerical solutions for specific gaseous systems. The spectro- 
scopic and correlation quantities needed for these calculations are available 
in [22,24]. 
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APPENDIX E 

EVALUATION OF CONSTANTS FOR STEADY LAMINAR FLOWS 

To determine the constants in Eq. (7.47), Eq. (7,49) is evaluated at 
5=0 and 5 = 1/4. To avoid excessive writing, the following notations are 
used (some of which are also used in Appendix C) 


5 ' = n , d 5 ' = dn 

b^. = 3 Uq^/2; c^. = 1 /b. = 2 / 3 Uq^ 


For 5 = 0, Eq. (7.49) reduces to 


aj + 1 +i (L/k) i H,.Uq. {/^[ai(Ti-2n3+n‘») 
2 i*l 0 


+ ag (r,2-2n3V)] A|(b.n)dT,} = 0 


(E.1) 


By defining u = b.n, Eq. (E.l) is expressed in an alternate form as 


ai{l+(L/k) I /^^[c^u-2(c.u)3+(c.u)‘t] ■Aj(u)du} 


+a2f(L/k) T H,. 
i=l 



(E.2) 


Now, by defining the following integral functions. 
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R 

R 

R 

R 


li 

2i 

3i 

4i 


Eq. (E.2) is expressed as 


0 

c 

I' 

0 

0 


u. 7^^ (u/) du. 

u? 7T: (u.) du. 
1 11 1 

u! 7T! (u.) du. 
111 ! 

U*t a: (u.) du. 
Ill 1 


(E.3a) 

(E.3b) 

(E.3c) 

(E.3d) 


ai aj + a2 «2 = ~1 


where 


ai = 1 t (L/k) Hi, (c,Ri, - 2 cf Rj, + R^,) 


a2 - (L/k) Hj, (c^Rj, - 2 c? R,, * c? R^,) 


For K = 1/4, Eq. (7.49) reduces to 


n. ^ .ii 

16 16 16 


^ ^ 1/4 1 

= i. (L/k) s {/ 0 (n) U>i (i - n)] dn 

2 i=l 0 ’ ^ 4 


fi/k "-)] du} 

^ ^ 4 


(E.4a) 

(E.4b) 

(E.4c) 


(E. 5 ) 
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By defining u = (i. - n) for the first integral and u = b. (n 

second integral, Eq. (E.5) is written as 


+ 11 aj + -1 a2 
16 16 16 


n b./4, 

= (L/k) E H,.{ f e(-- c.u) A!(u) du 
i=l 0 4 ^ ^ 

3b /4 , 

- / 0(i + c.u) 7T: (u) du 

0 4 '' '' 


By denoting d = c^u, the following relations are obtained from 


0(1- d) * a, (— - 11 d - -1 d2 + d3 + d‘») 

4 256 16 8 

+ cu(-i- - 11 d - — d2 - d3 + d*») 
256 16 8 

and 

0(1 + d) = a,(— + 11 d - -1 d2 - d3 + d**) 
4 256 16 8 

+ a« (- 1 - + _1 d - -1 d2 - d3 + d't ) 
256 16 8 


A combination of Eqs. (E.6) and (E.7) results in 


1 n C7 b./4 3b. /4 

ai(l E H J' (u) du - / ’ A'(u) du] 

K i=l 256 0 1 •’o 1 

11 l^i/i 3b. /4 

- — c.[/ uA'(u) du + / uA.'(u) du] 

16 1 0 1 0 1 


1) for the 
4 

(E.6) 

. (7.47): 

(E.7a) 

(E.7b) 
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9 b ./4 3 b . /4 _ 

- — C? [/ u2X'(u) du - / u2A'(u) du] 

8 ' 0 ' 0 1 

b ./4 3 b . /4 

+ c! [ / ^ u37T!(u) du + / ’ du] 

T o 1 0 1 

b ./4 3 b . /4 ,, 

+ cHC / ^ u^TTKu) du - / ^ u'fTT.'Cu) du]}- i±) 

I o'' 0 1 jLg 


b ,/4 


3 b . /4 


I II q ju./T 

+ 32 (i r H {- 1 - [ / ' 7 T[(u) du - / XHu) du] 

K i=l 256 0 1 ■’o 1 

o b ./4 3 b . /4 

- — c. [ / lu 7^1 (u) du + / u 7 ri (u) du] 

16 1 0 ■> 0 1 

1 b ^/4 3 b . /4 

- i c 2 [ / ^ u 2 X: (u) du - / ^ u 27 T:(u) du] 

8 ■> 0 ^ 0 1 

b ./4 3 b . /4 

+ c? [ / u 3 'ff^(u) du + / u 37 r^(u) du] 


b ./4 3 b. 74 

+ cV[ / ^ u'^7f! {u) du - / ^ u‘*7V!{u) 

I 0 1 0 ■> 


16 


11/16 


By noting that for any continuous function F(x) 


, 3 /** , ^ A/** . ^ , 3/4 , , 

/ F(x) dx = / F(x) dx + / F(x) dx 

0 0 1/4 


(E. 8 ) 


and defining 
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3b, /4 _ 

''{(“I) <'“1 (E-9a) 

^ “i ’^;<“l) "“i 

3b, /4 

^3i -I, “l ^ (E.9C) 

3b, /4 

^41 ' ! ,. “1 ‘‘“I (E-M) 

D^./4 
"i'''* 3 

^5i ' “i ‘'“l (E-9b) 

3b, /4 

^6l ' “1 ^<“f) "“1 (E-9f) 

3b /4 ^ 

Syi ' / “i ’fi('-i) <*“1 (E-99) 

Eq. (E.8) can be written as 

3l 03 + 32 04 ~ - 11/16 (E.lOa) 


where 
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n 

03 = (L/k) [(57/256) + (11/16) (S2i+S3^) 

- (9/8) c2. - c? (Sg^+Sg.) + cl S7/] + 11/16 (E.lOb) 

n 

04 = (L/k) (9/256) $ 1 ^ + (3/16) c,- (S2^+S3^) 

- (1/8) c? S4. - c? (Sg^+Sg.) + e‘tS7.] + 3/16 (E.lOc) 

By solving Eqs. (E.4a) and (E.lOa) simultaneously, there is obtained 
the results for constants a^ and aj as 

ai = (11 «2 - 16 04) /DEN (E.lla) 

ag = (16 03 - 11 oi)/DEN (E.llb) 

where 

DEN = 16 (0204 “ 0203) (E.llc) 

The governing equation for the large path length limit is Eq. (7.44a) 
for which the solution is also given by Eq. (7.47). For the large path 
length limit, Eq. (7.47) is expressed in the form of Eq. (7,52) which is 
represented here as 

e(s) = bi (C - 2?3 + cM + b2 (5^ - 2«V+ S'*) (E.12) 

Thus, 

e'(5) - bid - 652 + 453) + bi (2| - 6f2 + 453) (E.13) 


c-^ 
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A substitution of Eqs. (E.12) and (E.13) in Eq. (7.44a) gives 

bj (1 - 6c2 + 453) + (25 - 652 + 453) - 2 (3g2 - 253) + 1 

= M f^0(5') d57(5-5') (E.14) 

0 

For 5 = 0, Eq. (E.14) reduces to 

b. + 1 = - M [9(5 ')/5'] d5' (E.15) 

0 

Upon substituting for 0(5') from Eq. (E.12) into Eq. (E.15), the integrals 
can be evaluated in closed form and there is obtained 

bj 01 + b2 &2 ~ ~ ^ 

where 

81 = 1 + (7/12) M; 82 = ^1/12) M 

For 5 = 1/4, Eq. (E.14) reduces to 

^b, +-^b, +^=M r%(5') dc7(-- ?■) (E.17) 

16 16 16 0 4 

By substituting for 0(5') in this equation, another relation between bj and 
b2 can be obtained in terms of 83 and 84. But, this appears to involve the 
evaluation of a few improper integrals. Thus, this approach is abandoned in 


(E.16a) 

(E.16b) 
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favor of an alternate procedure discussed below. 

The solution in the large path length limit can be obtained from the 
general solution by evaluating the integrals R. and S. in the limit of large 
path length. Since in this limit A'(u) = 1/u, the integrals in Eqs. (E.3) 
and (E.9) are evaluated to obtain 

hi = ‘^i’ ‘^2i = “^Si “ *^4i " 

= in (3); $2^ = b./4; S3. = 3b./4; S^. = b\/^ ; 

S5. = b3/192; Sg. = 9b? /64; S7^ = Sb'jM (E.l8b) 

From Eqs. (E.4) and (E.18a), there is obtained for the large path length 
limit 


bi 61 + b2 02 ■ ~ 1 


(E.19a) 


Bl = 1 + (7/12) M; 02 = (1/12) M (E.19b) 

which is the same result as given by Eq. (E.16). From Eqs, (E.IO) and 
(E.18b), one obtains in the large path length limit 


bi 03 + b2 04 ” " 11/16 


(E.20a) 


where 


03 - 11/16 + M [(57/256) jtn(3) + 65/192] 
= 11/16 + 0.583 154 559 M 


(E.20b) 
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3 ^ = 3/16 + M [(9/256) £n(3) + 17/192)] 
= 3/16 + 0.127 164 755 M 

The solution of Eqs. (E.19a) and (E.20a) yields 


bi = (1132 - 16 34 ) /BOTTOM 


b 2 = (1633 - 11 6i)/B0TT0M 


where 


BOTTOM = 16 ( 3 i 34 - 6263 ) 


With bi and 62 known, the solution for the temperature distribution 
obtained from Eq. (E.12). 


(E.20c) 

(E.21a) 

(E.21b) 

(E.21c) 

is 
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APPENDIX F 

INTEGRAL FUNCTIONS FOR TRANSIENT LAMINAR FLOWS 

For convenience and use in the computational procedure, the following 
definitions are employed in expressing the relations for the integral 
functions: 


b. = 3 Ujj/2; c^ = 1/b.; C(5) = l/(?-?2) 

X = ? - 253 + s'*; y = ?2 _2g3 + ^>*; z = 12 (s-s2) 

a(0 ^ 1 - 6?2 + 4s 3; b(s) = 2 ? - 652 + 4s3 


Various integrals are defined as follows: 


2 _ 


Ji(S) = 3 N + C(S) (2 +1 (L/k) Z Hi. uL{/ (?•-?' ) AV [b.(S-S-)] d?' 

4 1=1 0 ' ' 


+ / U'-S' ) [b (S'-S)]d5'}) 
S 


(F.l) 


J2U) = CU) (2 +- (L/k) I Hi.Uq.{/ (l- 2 g') A! [b.(g-5')]dc‘ 

2 i “1 ft 


/ (l-2c') X' [b.(s'- 5 )]dc'}) 
S 


(F. 2 ) 


q " 5 

JaU) = 2 + 3 Nx +_ (L/k) E H -u2 // x(f ') X' [b (g-5')] d^' 
^ ^_1 n 01 ‘■' 0 T ^ 
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1 

+ / xU')A? [b.(€'-?)]d 5 '} (F.3) 

q n c 

04 ( 5 ) = (z-2) + 3Ny + 1 (L/k) z H u2 {/ y(g') A!! [b.C?-? ' )]d? ' 

4 * 1 =]^ 0 

1 

+ y{5') [b.(5-5 ')]d€'} (F.4) 

on 5 _ 

J 5 (C) = 2 +- (L/k) z H, u .{/ a(?V) A1 [b dC 

2 1=1 0 • ' 

1 

- a(? ' ) A^' [b-j (? ■ < ) ]d5 ' } ( f .5) 

/ on C 

JeU) = (z - 2) +- (L/k) E H u {/ b(r') A! Cb.(c-?')] dc* 

2 i=l 01 0 ^ 1 

1 

- / t)(6') a; [b,(5'-e)]d5'} (F.6) 

1 

07 ( 5 ) ^ 3N + C( 5 ) [2 + M / (1-2C) d?7(?-C')] 

0 

= 3N + C(c) f2 + M [2 + (2f-l) in (ill)]} (F.7) 

C 

JsC?) = <> 7 ( 5 ) - 3N (F.8) 



